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Abstract 
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CO ' We present the complete four-loop four-point amplitude in N = 4 super- Yang-Mills theory, for a 

00 

■ general gauge group and general D-dimensional covariant kinematics, and including all non-planar 

o : 

contributions. We use the method of maximal cuts — an efficient application of the unitarity 
method — to construct the result in terms of 50 four-loop integrals. We give graphical rules, valid 
5_i ■ in .D-dimensions, for obtaining various non-planar contributions from previously-determined terms. 

We examine the ultraviolet behavior of the amplitude near D = 11/2. The non-planar terms are 
as well-behaved in the ultraviolet as the planar terms. However, in the color decomposition of the 
three- and four-loop amplitude for an SU(N C ) gauge group, the coefficients of the double-trace 
terms are better behaved in the ultraviolet than are the single-trace terms. The results from this 
paper were an important step toward obtaining the corresponding amplitude in N = 8 supergravity, 
which confirmed the existence of cancellations beyond those needed for ultraviolet finiteness at four 
loops in four dimensions. Evaluation of the loop integrals near D = 4 would permit tests of recent 
conjectures and results concerning the infrared behavior of four-dimensional massless gauge theory. 
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I. INTRODUCTION 



In recent years, scattering amplitudes have become an important tool for studying fun- 
damental issues in gauge and gravity theories. For example, for the maximally supersym- 
metric four-dimensional super- Yang-Mills theory (TV = 4 sYM) in the planar limit, we have 
recently seen an all-order resummation ansatz proposed by Smirnov and two of the authors 



(BDS) , strong-coupling results from Alday and Maldacena [3| and the identification 
of of a new symmetry — dual (super) conformal invariance j^f]]. Together, these results 
have opened a new venue for studying the holographic AdS/CFT correspondence. Further- 
more, they suggest the remarkable possibility that planar amplitudes of M = 4 sYM may 
ultimately be determined exactly [7j. In the maximally supersymmetric four- dimensional 
supergravity theory (J\f = 8) jg], the remarkable ultraviolet behavior of multi-loop graviton 
scattering amplitudes (q- 12] has challenged the widely-held belief that it is impossible to 
construct a perturbatively consistent point-like theory of quantum gravity. 

In this paper we compute and analyze the four-loop four-point amplitude in M = 4 sYM, 
for an arbitrary non-abelian gauge group G. The amplitude can be graphically organized 
into planar and non-planar contributions. For the case of a special unitary gauge group, 
G = SU(N C ), the planar contributions dominate the limit in which the number of colors, N< .,, 



13, 



13|- 



tends to infinity. The planar terms are much simpler, and were computed previously 
The non-planar contributions are the subject of this paper. 

Scattering amplitudes in gauge theory and gravity have a much richer structure than 
that revealed by Feynman diagrams. One striking example is Witten's observation that tree 
amplitudes are supported on curves in twistor space [ljj]. At the multi-loop level, such struc- 
tures have first been identified in theories with maximal supersymmetry. In M = 4 sYM, 
the first hint of powerful relations between higher- and lower-loop amplitudes was the "rung 
rule" relation 15|, ]l6] for constructing particular contributions to four-point amplitudes. 
More remarkably, an iterative structure was uncovered in the dimensionally-regularized pla- 
nar two- loop four-gluon amplitude [lj], using the explicit values of the loop integrals in an 
expansion around D = 4 17j. This iterative structure holds at three loops as well; it led to 
the all-loop BDS ansatz [2 ] for maximally-helicity- violating (MHV) amplitudes. For four ex- 
ternal states, the ansatz is almost certainly correct, given its verification at strong coupling 
by Alday and Maldacena js| using the AdS/CFT correspondence. For four or five external 
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legs, the assumption of dual conformal invariance, together with the structure of the infrared 
divergences, completely fixes the amplitudes' dependence on the kinematics {(J, making it 
very likely that the BDS ansatz is correct. In addition, a variety of explicit four- and five- 
point computations confirm it through three loops in dimensional regular izat ion 
and with a Hi ggs reg ulator j^oj ]. However, for six or more external states, the BDS proposal 



is incomplete 2lN24|. To clarify the structure of the additional "remainder" terms, it will 



undoubtedly be important to carry out further computations at both weak and strong cou- 
pling. Positive steps have been taken in this direction recently [7, [25!^ 27] . 

Much less is known about non-planar contributions to M = 4 sYM amplitudes. As a 
step in this direction, the principal aim of this paper is to construct the complete four-loop 
four-point amplitude in M = 4 sYM, including non-planar contributions, using the unitarity 
method 28( and various refinements of it. We express the result in terms of 50 distinct loop 
momentum integrals, each with nontrivial numerators. (One of the 50 gives a contribution 
to the amplitude that vanishes after integration.) At three loops, the analogous expression 



requires only nine distinct integrals 10 



29]. 



Using these representations of the three- and four-loop amplitudes, we can study their 
ultraviolet (UV) properties, in particular the critical dimension D C (L) in which the first UV 
divergence appears, as a function of the loop number L. Based on information from some 
of the unitarity cuts, a formula for D C (L) was suggested in ref. 16J (see eq. (15. ip below). 
This formula corresponds to counterterms in higher dimensions of the schematic form V 2 F 4 , 
where F is the Yang-Mills field strength, T> is a covariant derivative, and the precise color 
structure is not yet specified. This form for the counterterms was later argued for by Howe 
and Stelle [30j, based on harmonic superspace [311 ] . 



Here we confirm this expected behavior for the so-called single-trace terms in the ampli- 
tude for G = SU(N C ), which correspond to counterterms of the form Ty(T> 2 F a ). However, 
we also find striking cancellations in the double-trace terms in the amplitude, starting at 
three loops, such that counterterms of the form Tt(T> 2 F 2 ) Tr(F 2 ), etc., are absent. These 
results were first reported in refs. 32J. These cancellations have been discussed using the 



pure spinor formalism in string theory 



331 ] and fie 



d theory 



12J, as well as from the point 



341 ] . Because of contamination from the 



of view of algebraic non-renormalization theorems 
closed-string sector, the string-based arguments apply only to the double-trace terms with 
the most factors of A c , namely A^ _1 ; the field-theory arguments [3] are more general in 
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this regard. Here we will rearrange the color structure of the amplitudes in order to make 
manifest the additional cancellations in the double-trace terms, including all powers in N c , 
at three and four loops. 

The present paper also provides the key input for the unitarity-based computation of the 
four- loop four-graviton amplitude 11] in J\f = 8 supergravity [8j . On general grounds, the 
unit arity cuts of gravity loop amplitudes can be obtained from those of Yang-Mills ampli- 



tudes 



161 ]. In a first step, generalized unitarity allows gravity loop amplitudes to be de- 



composed into gravity tree amplitudes. Then the tree-level Kawai, Lewellen and Tye (KLT) 



relations 



35] can be used to express the gravity tree amplitudes in terms of gauge-theory tree 



amplitudes. Alternatively, one can use the new diagrammatic numerator relations between 



gravity and gauge-theory tree amplitudes developed by three of the authors [36]. Thus the 
gravity cuts are expressed as bilinear combinations of gauge-theory cuts. All the information 
needed for the gravity computation may be found by cutting the gauge-theory amplitude. 
In the case of M = 8 supergravity, all the sums over supersymmetric particles running in the 
loops are automatically performed in the course of the M = 4 sYM computation. Because 
gravity does not involve color, the non-planar contributions are just as important as the 
planar ones, and the full non-planar M = 4 sYM amplitude is required. 

Multi-loop M = 8 supergravity amplitudes have revealed UV cancellations not anticipated 
from earlier superspace power-counting arguments. The explicit UV behavior of the three- 
and four-loop four-point amplitudes 10|, 111] exhibits strong cancellations, beyond those 
needed for UV finiteness in four dimensions. Moreover, an analysis of certain unitarit 
cuts demonstrates the existence of novel higher-loop cancellations to all loop orders 
based on the one- loop "no-triangle property" [37H42]. These results support the proposal 



ity 

I 



that M = 8 supergravity might be a perturbatively UV finite theory of quantum gravity. 
String dualities [43] and non-renormalization theorems 44J have also been used to argue 
for both UV finiteness of M = 8 supergravity, and for a delay in the onset of divergences, 
although difficulties with decoupling towers of massive states 45] and technical issues with 
the pure-spinor formalism 12j, |46|, |47] may affect these conclusions. 

Another important reason to study the non-planar contributions to M = 4 sYM am- 
plitudes is to investigate whether the resummation of the four-point amplitude to all loop 
orders [l ^] can be accomplished in some form for non-planar amplitudes as well. However, 
given the more complicated color structure of non-planar contributions, and the probable 
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lack of integrability for sub leading terms in the 1/N C expansion, it is not clear whether such 
a generalization can exist. What is clear is that the infrared- singular behavior would have 
to be understood first, before the behavior of infrared-finite terms could be addressed. At 
sub leading orders in 1/N C , the soft anomalous dimension matrix T$ controls infrared sin- 



48 



491 ] . The matrix structure becomes trivial in the 



gularities due to soft-gluon exchange 
planar, or large- N c , limit j^j. 

Surprisingly, the two-loop soft anomalous dimension matrix in massless gauge theory 
is proportional to the one-loop one 50|. The proportionality constant is determined by 



the cusp anomalous dimension 



491 ] . This result was later understood to be a consequence 



of an anomalous symmetry of Wilson-line expectation values under the rescaling of their 



velocities 



51 



52| . It has been conjectured 



29 



51 



52J that the proportionality might persist 



to all loop orders. However, velocity rescaling a 



51 



one, even combined with other constraints, 



53| , is not powerful enough to determine the 



such as collinear factorization of amplitudes 

form of Ts beyond two loops j^], except for the matter-dependent part at three loops [55]. 
In principle, the soft anomalous dimension matrix for M = 4 sYM can be determined at 
three and four loops, for the case of four external massless states, using the results in ref. 29] 
and in this paper. At three loops, this computation would also determine Ts for any massless 
gauge theory, because the matter-dependent part is known [55]. 

To perform such a determination, the dimensionally-regularized loop integrals entering 
the amplitudes, for D = 4 — 2e, would have to be expanded in a Laurent series around 
e = 0. If this can be accomplished, then it is relatively straightforward to compare the 
results with fixed-order formulae 50, [56] in order to extract Ts- The Laurent series for the 
L-loop amplitude begins at order l/e 2L , and the L-loop soft anomalous dimension matrix 
first enters at order 1/e in the expansion. 

At order 1/e 2 , the L-loop cusp anomalous dimension 7^ appears. While this quantity 
is known in the planar limit of M = 4 sYM through at least four loops [3, , and very 
likely to all loop orders 58[, the sub leading-color terms are not known beyond three loops. 
If they are nontrivial at four loops, it would indicate the violation of "quadratic Casimir 



scaling", or 7^ oc for G = SU(N C ), which holds through three loops 



Such 



a violation is allowed by group theory, beginning at four loops, and hinted at by strong- 



coupling considerations 60|. However, it would be very useful to compute the subleading- 
color terms in 7^ explicitly, as it has been conjectured that quadratic Casimir scaling will 
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continue to hold beyond three loops 51 |. 

Unfortunately, both these tests will have to await improved techniques for the evaluation 
of the Laurent expansion in e of three- and higher-loop non-planar integrals. The non-planar 
integrals required for three-loop form factors (one off-shell leg and two on-shell massless ones) 



have been evaluated 



611 ]. but not yet those for on-shell scattering of four massless states at 



three loops (let alone four loops). 

In order to determine the four-loop amplitude in terms of a set of loop integrals, we use 
the unitarit y m ethod 281 ]. in particular the information provided by generalized unitarity 



cuts 



53 



62j-[64|. The unitarity method takes advantage of the fact that tree-level amplitudes 



are much simpler than individual Feynman diagrams. One can exploit structures that are 
obscure in Feynman diagrams but visible in on-shell amplitudes. 

For planar M = 4 sYM, the observation that the loop amplitudes can be expressed in 



terms of a restricted set of "pseudoconformal" integrals [4j, [13j, |64j, [65| helps streamline the 
construction of a candidate expression, or ansatz, that is consistent with the unitarity cuts. 
For a planar graph, dual coordinates Xi can be associated with nodes of the dual graph; 
differences of pairs of Xi are identified with momenta of the internal or external lines of the 
original graph. Loosely speaking, pseudoconformal integrals are integrals that are invariant 
under conformal transformations of the dual coordinates. It has not yet been proven that all 
planar loop amplitudes in M = 4 sYM are expressible in terms of pseudoconformal integrals 
(especially for non-MHV helicity configurations for six or more external states). However, 
in any given loop amplitude we can directly confirm that no other integrals appear, once 
the pseudoconformal ansatz satisfies all the unitarity cuts. This observation helped simplify 
recent calculations of the planar five-loop four-point 64j and two-loop six-point 23| MHV 
amplitudes in M = 4 sYM. 

In contrast, for non-planar terms in the amplitude, there is currently no useful definition 
of dual conformal symmetry. In addition, non-planar contributions are inherently more 
intricate than planar ones. Nevertheless, planar and non-planar contributions do appear 
to be linked via identities between diagram numerators, whose structure is similar to the 



Jacobi identity obeyed by color factors 



36 



66 



671 ] . Very recently, this structure has been 

. Indeed, for 



observed to hold at the three-loop level, with no cut conditions imposed 
the three-loop four-point amplitude in M = 4 sYM, one of the planar diagrams is sufficient 
for determining all the non-planar ones. It would be interesting to study these properties 
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further, in particular to determine the restrictions they place on the four-loop amplitude 
presented here. We leave this interesting question to the future. 

A powerful way to determine large classes of both planar and non-planar contributions 



is through graphical rules that capture some o 



The first of these rules is the rung rule 15|, [l6|. The Jacobi-like relation between planar 



the features of generalized unitarity cuts. 



and non-planar contributions can also be implemented as a graphical manipulation 
There is also a "box-substitution rule" J64J, which we will generalize here. Some related 
graphical identities may be found in ref. [69|. The rules given here reduce the determination 
of contributions containing four-point subdiagrams to a few simple manipulations. 

The graphical rules presented in this paper capture many contributions, but not al 
them. In order to determine the missing ones, we used the method of maximal cuts 



29 



64l | . For four-point amplitudes in D > 4, maximal cuts are those generalized unitarity 



cuts in which only three-point tree amplitudes appear, connected by cut propagators, i.e., 
they contain the maximum number of cut propagators. One first constructs a candidate 
expression (ansatz) for the loop amplitude that is equal to each maximal cut, when evaluated 
for the appropriate cut kinematics. We refer to this procedure as "matching" the ansatz to 
the maximal cuts. Next one adjusts the ansatz to match as well the next-to-maximal cuts 
in which a cut propagator is removed from between two of the three-point trees to form a 
four-point tree. This captures all contributions with a single four-point vertex, or "contact 
term" . Then one matches to further near-maximal cuts with more canceled propagators. 

In order to construct the complete four- loop amplitude, we used a large number of max- 
imal and near-maximal cuts. To ensure that no contributions were dropped, we then evalu- 
ated a set of 13 "basis cuts" (not counting permutations of legs) which suffice to determine 
any massless four-loop four-point amplitude. We verified that our answer matches all of 
these cuts. Many of the cuts were evaluated using four-dimensional intermediate momenta, 
so that powerful helicity and supersymmetry methods could be exploited. This leaves open 
the possibility that some terms could be missed, which vanish when the cut momenta are 
restricted to four dimensions. For this reason, we performed a large number of consistency 
checks, as described in section III Bl 

The unitarity method also allows us to use an on-shell superspace formalism to sum 
over states crossing cuts. On-shell superspaces involve only physical states and, for our 
purpose, they are generally far simpler than their off-shell cousins. A number of years ago 
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Nair presented an on-shell superspace 70J for MHV tree amplitudes in M = 4 sYM. More 
recently, this superspace has been extended to any helicity and particle configuration. For 



the computations of this paper we follow the MHV generating-function approach 7lH73| . as 



organized in ref. [iJ] for use in multi-loop calculations. We also exploit specific super-sum 



results for cuts involving next-to-MHV tree amplitudes 73J. A related procedure 



42 



75 



76| for covering any helicity and particle content employs the momentum shifts used by 



Britto, Cachazo, Feng and Witten to derive on-shell recursion relations 



771 ]. extended to 



shifts of anti-commuting parameters. One-loop examples that use an on-shell superspace in 



conjunction with the unitarity method may be found in refs. 
examples, including four-loop ones, are given in refs. 



73 



74|. 



42 



72 



761 ] . Various higher-loop 



This paper is organized as follows. In section [TT] we describe the general structure of 
multi-loop amplitudes. We also recall the specific form of four-point amplitudes in M = 
4 sYM from one to three loops, as well as the planar terms at four loops. We give a brief 
overview of the techniques used to determine the amplitudes. In section III II we describe 
in more detail various tools for determining the non-planar contributions. In section IIVI 
we present the complete four-loop four-point amplitude in terms of a set of 50 integrals. 
The ultraviolet divergence properties of the four-point amplitude through four loops, in the 
critical dimension D C (L), are discussed in section [V] In section IVIl we compute the leading 
UV divergence for the double-trace terms at three loops, which appears at D = 20/3 (in 
contrast to the single-trace terms, which first diverge in D c (3) = 6). In section IVIIl we give 
our conclusions and prospects for the future. In appendix we present a sample evaluation 
of a nontrivial non-planar cut. In appendix [B] we provide various representations of the 
color factors appearing in the amplitudes. In appendix [C] we collect the numerator and 
color factors for the 50 integrals entering the four-loop amplitude. 



II. STRUCTURE OF MULTI-LOOP AMPLITUDES 



Loop amplitudes in M = 4 sYM exhibit remarkable simplicity for a gauge theory. Using 
the unitarity metho d |28| . a large variety of amplitudes have been constructed through five 



loops [13 



64 . l78| in terms of loop-momentum integrals. Indeed, the struc- 



ture of the planar four-point amplitude is simple enough that an all-loop order resummation 
is possible [2l, |3[. This simplicity in the planar sector has been understood in terms of a new 
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symmetry dubbed "dual conformal symmetry" [4j, |5|, [65(, which is intimately connected to 



integrability 



an 



791 ] . In this paper we focus on the non-planar contributions to M = 4 sYM 



amplitudes. Although they are much more intricate and less well understood than the 
planar amplitudes, their structure is still remarkably simple, especially when compared to 
amplitudes in theories with fewer supersymmetries. 

In this section we begin by describing the color and parent-graph organization of multi- 
loop amplitudes, including a review of the results for the lower-loop and planar four-loop 
four-point amplitudes in M = 4 sYM. Then we turn to a brief review of the unitarity method. 



A. Color and parent-graph decomposition 

For gauge group G = SU(N C ), the leading-color (planar) terms are particularly simple. 
They have essentially the same color structure as the corresponding tree amplitudes. The 
leading-in-A c contribution to the L-loop n-point amplitude may be written as, 

L 

TrCTMu . . . TM"))A( L )(a(l), a(2), . . . , a(n)) , (2.1) 

n 

where the T ai are generators in the fundamental representation of SU(N C ), with adjoint 
color indices a,, and the sum runs over non-cyclic permutations of the external legs. In this 
expression we have suppressed the (all-outgoing) momenta ki, as well as polarizations and 
particle types, leaving only the integer index i as a collective label. This decomposition holds 
for all particles in the gauge super-multiplet, as they are all in the adjoint representation. 
The color-ordered (or color-stripped) partial amplitudes carry no color indices; they 
depend only on the kinematics, polarizations and particle type. At leading order in the 1/A C 
expansion they can be expressed solely in terms of planar loop integrals. 

For the complete amplitude for a general gauge group G, including all non-planar contri- 
butions, the parent-graph decomposition, 

A (L) = g 2L + n-2 J- OiGIt, (2.2) 
i G parent 

is more convenient than the color-trace representation. The parent graphs are cubic graphs 
- graphs containing only three-point vertices. Momentum is conserved at each vertex. 
Every graph specifies simultaneously a combinatorial factor <2j, a color dressing Cj and a 
Feynman loop integral /j. 

9 



. Planar 
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(4vr) 2 - 



FIG. 1: Rules for obtaining the color factors associated with cubic parent graphs. The roman 
letters a, b, c are color indices. The sign of each structure constant f abc is fixed by the clockwise 
ordering of each vertex, as drawn. 

The Ci are written in terms of group structure constants. The contractions of group 
indices are encoded by the graph i using the rules given in fig. [TJ More precisely, the color 
factors Ci are obtained by dressing each three-vertex of the parent graph with structure 
constants, normalized as, 



where / are the standard structure constants of the gauge group G, and the hermitian 
generators T a are normalized via Tr (T a T b ) = 5 ab . The f abc should follow the clockwise 
ordering of the legs for each vertex of the parent graph. This clockwise ordering is important. 
For example, if one redraws a graph so that an odd number of three-point vertices have their 
ordering reversed, then the signs of both C» and l{ should be flipped in tandem. 

Each 'parent integral for the L-loop four-point amplitude is a Feynman integral with the 
following general structure and normalization, 



where k m , m — 1, 2, 3, are the three independent external momenta, £j are the L independent 
loop momenta, and l n are the momenta of the (3L + 1) propagators (internal lines of the 
graph i), which are linear combinations of the £j and the k m . As usual, d D £j is the D- 
dimensional measure for the j th loop momentum. The numerator polynomial Ni(£j,k m ) is 
a polynomial in both internal and external momenta. 

Unlike the decomposition using color traces, the parent-graph decomposition is not 
unique, due to contact terms. Contact terms are contributions to the amplitude that lack 



abc 



Tr([T a , T b ]T c ) , 



(2.3) 




(2.4) 
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one or more propagators, relative to the parent graphs. If the contact terms were allowed 
to contribute as isolated integrals, they would correspond to graphs containing quartic or 
higher-order vertices. Here we will absorb all contact terms into parent integrals, by mul- 
tiplying and dividing by the missing propagator or propagators, so that the corresponding 
term in Ni(£j, k m ) will contain factors of 1^, which we refer to as inverse propagators. How- 
ever, because the associated color factors Ci can be expressed as linear combinations of other 
color factors via the Jacobi identities, there is an ambiguity in choosing the specific parent 
graph to absorb a given contact term. Although there are many valid choices, particular 
choices can reveal nontrivial structures or symmetries. 

For M = 4 sYM amplitudes, the freedom in assigning contact terms to parent graphs 
can be exploited to remove all graphs with nontrivial two- or three-point subgraphs, as was 
done at three loops 29|. While supersymmetry and gauge invariance may be used to show 



the all-order consistency of this condition on the parent graphs, a more direct approach was 



taken in ref. 



29J, by verifying that an ansatz in this class is compatible with all unitarity 



cuts. Here we will take the same approach at four loops. 

2. . 3 



FIG. 2: The box diagram is the only parent graph at one loop (up to permutations of the external 
legs). 



At one loop, the structure of the M = 4 sYM four-point amplitude is especially simple. 
We modify eq. (I2.2p slightly by extracting an overall prefactor, and write the result as, 



A 



(i) 



box rbox ; 
1234 1 I 



Sl2,S 2 3j 



(2.5) 



S4 



where g is the gauge coupling. The prefactor K is defined by 



K = /C(l, 2, 3, 4) = s 12 s 23 A?*(l, 2, 3, 4) 



(2.6) 



where Si 2 = (k\ + k 2 ) 2 , s 23 = (k 2 + k-j) 2 , and the ki are the external momenta. It contains 
all information about the four external states. 
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°1234 



The unique parent graph at one loop is the box diagram shown in fig. [2j For external 
legs ordered 1234, the box color factor is 

j-aiaea?, ja^aiaa j-asagar ja^a^as ^ <j\ 

where we sum over repeated indices. This form is valid for any gauge group. Finally, 
J box (si2, S23) is the one-loop box integral, 

(2.8) 



/ b ° X (Sl2,S 23 ) 



{27r) D p 2 (p - ktfip - h - k 2 ) 2 (p + h) 2 ' 



The sum in eq. (12. 5ft runs over the 24 permutations of external legs {1, 2, 3, 4}, denoted by 
S4. The permutations act on both the momentum and color labels. The prefactor of 1/8 
accounts for an eightfold overcount in the permutation sum, which we leave in to make it 
slightly easier to generalize to higher loops. 

In eq. fl 2 . 6 j) . A^ ee (l, 2, 3, 4) stands for any M = 4 sYM tree amplitude in the canonical 
color order. In four dimensions, a compact form of this object can be written down using 



anti-commuting parameters rj in the on-shell superspace formalism 42 

i<5 (8) (A!?7i + \ 2 r]2 + A3773 + A4774) 



70 



73 



75 



Af e (l,2,3,4) 
£(1,2,3,4)1 



D=A 



(12) (2 3) (3 4) (41) 
[12] [3 4] 



<5 (8) (Ai?7i + \ 2 r] 2 + A37/3 + A4774) 



(2.9) 
(2.10) 



lz? = 4 (12) (3 4) 

It is not difficult to verify that ([12] [3 4])/((l 2) (3 4)) is symmetric under exchange of any 
two legs, and that JC is local. Related to this, JC represents the color-stripped four-point 
(linearized) matrix elements of the local operator TrF 4 , plus its supersymmetric partners. 
Therefore K, is a natural prefactor to extract from the four-point amplitude in M = 4 sYM. 



S12 




/ 



(P) 



(NP) 



V 



FIG. 3: The planar and non-planar parent graphs contributing to the two-loop four-point ampli- 
tude. The prefactors s\2 are the integral numerators and N^ p \ respectively. 

At two loops, the full M = 4 sYM amplitude is given by a similar permutation sum as 
for the one- loop case (I2.5P flol . Q, 
1(2) 



at = -/^Y, p24 /(p) (wad + C^ (NP) (%,%) 

Si 



(2.11) 



12 



The planar and non-planar double-box integrals, displayed in fig. [31 are defined by 

d D pd D q 



I {P \s 12 ,s 23 ) = {-if 



Sl2 



(27T)2^ p2 (p _ fcl )2 (p _ kl2) 2 (p + g) 2 g 2 (g _ fc4 )2 (g _ ^2 



d D pd D q 

(27r) 2jD p 2 (p - k 4 ) 2 (p + g) 2 (p + q + h) 2 q 2 (q - ki) 2 (q - k 12 ) 2 ' 

(2.12) 



with ki 2 = k\ + k 2 and = k% + £4. The permutation sum again runs over S4 and acts on 
both momentum and color labels. Because both graphs in fig. [3] have a fourfold symmetry, 
the permutation sum overcounts by a factor of four. As before, we include this overcount 
and divide by an overall symmetry factor. In a form valid for any gauge group, the color 
factors of the planar and non-planar graphs, with legs ordered 1234, are, 



^y(P) _ ja\CL6CL5 jaia~ia§ ja^agas jaiaioag jajasan ^aganaw 
£,(NP) _ jaia 6 a 5 ja 2 a 7 a 6 Ja 3 aga s ^04011010 Ja 7 aga 1Q Ja 5 ayan 



(2.13) 



At three loops, the fully color-dressed three-loop four-point Af = 4 sYM amplitude is 



given by 



10 



22| 



A {3) 



4 Si 

+ 2C (c) / (e) + 2C (f) / (f) + 4C (g) / (s) + |C (h) J (h) + 2<7 (i) / (i) 



• (2-14) 



In this case, the integrals I^(si 2 , S23) are D-dimensional loop integrals corresponding to the 
nine graphs shown in fig. HJ using eq. (12. 4p with the numerator polynomials iVj displayed 
next to the diagrams. In fig. [4] the (outgoing) momenta of the external legs are denoted by 
ki with % = 1,2,3,4, while the momenta of the internal legs are denoted by U with i > 4. 
For convenience we use the following shorthand notation, 



(ki 4~ kj) 2 
(ki + lrf 
(h + lj) 2 



i,j < 4 
i < 4 < j 
4 < i,j 



Ti 



(2.15) 



In the three-loop case, some of the numerator polynomials contain squares of loop momenta, 
which could be used to collapse propagators and generate contact terms. (The three-loop 
four-graviton amplitude in M = 8 supergravity can be written 29| in a form very similar 
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FIG. 4: The nine parent graphs for the three-loop four-point AA = 4 sYM amplitude. The pref- 
actor of each diagram is the numerator polynomial iV^'. The kinematic invariants are defined in 
eq. (j2T5|) . 



to eq. f)2.14p . except that there are no color factors, and the numerator factors in the loop 
integrals are of course different.) 

The color factor associated with each integral in the three-loop amplitude is easy to 
write down from the parent graph, following fig. [TJ The expression ( I2.14p is valid for any 
gauge group G and any dimension D, as verified by a direct evaluation of the color-dressed 
cuts [29]. 

In section HVl we will present the complete four-loop four-point amplitude of M = 4 sYM, 
using the same type of parent-graph decomposition. As we will demonstrate, this amplitude 
can be decomposed into 50 distinct parent integrals, corresponding to cubic graphs with 
no nontrivial two- or three-point subgraphs. As explained earlier, contact terms are incor- 
porated as numerator factors containing inverse propagators. Thus we write the four-loop 



14 



'12 



2 

<Sl2<S56 



(d) 



Sl2-S45'536 



i 1 










(b) 


\ 
y 








5 


1 r 


(e) 



2 

s 12 s 45 











5 





•Sl2-S45'S26 



(c) 



1 (f) 



FIG. 5: Planar contributions to the four-loop four-point amplitude. 
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FIG. 6: Contact terms in the planar four-point amplitude, which can be absorbed into the numer- 
ator of graph (f) of fig. [5j (Note that here (f2) and (d2) refers to the labeling used in ref. [ljj; 



unlike ref. 



131 ] we choose to write out explicitly the relative factors of these integrals with respect 



to diagram (f).) 



50 



four-point amplitude as, 

S 4 i=l 

where the usual prefactor JC is common to all terms, the are combinatoric symmetry 
factors and the color factors. The integrals Ii are specified by the propagators associated 
with the parent graph, and by the numerator polynomials iVj. Each numerator polynomial 
is subject to various constraints. After accounting for the four powers of external momenta 
in K, in eq. (12.161) . dimensional analysis implies that the numerator polynomial Nj (£j , kj) is 
of degree 6 in the momenta. Moreover, consistency with the known power-counting [lq . l30l ] 
requires that the numerator polynomials for M = 4 sYM have a maximum degree of 4 in 
the loop momenta. 
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The planar contributions to the four-loop four-point amplitude were presented earlier 13| , 
although they were given in the color-trace decomposition, rather than the graphical de- 
composition used here. Out of the eight different types of planar integrals present in the 
amplitude, the six parent graphs are shown in fig. (The combinatoric factors are not 
included in the figure.) The two additional contributions are contact terms and are shown 
in fig. [6j A convenient choice for absorbing these contact terms is to assign both of them to 
the parent graph (f ) of fig. [5j by incorporating inverse propagators in their numerators. The 
result, after combining two different permutations of diagram (f), can be found in fig. 
diagram (28). 



B. Unitarity method 



The unitarity method provides an efficient framework for systematically constructing and 



verifying the expression for any massless multi-loop amplitude. This methoc 



various refinements, has already been described in some detail elsewhere 
64, 



16 



along wit 



28 



53, 



62- 



69|. Here we summarize those points directly salient to our construction of multi-loop 



amplitudes in TV = 4 sYM. 




FIG. 7: Examples of maximal and near-maximal cuts used to construct a complete ansatz for the 
four-loop N = 4 sYM amplitude. Cuts (a) and (b) have the maximal 13 cut conditions. Near- 
maximal cuts (c) have 12 conditions, and (d) and (e) have 11 conditions. No cuts with fewer than 
11 on-shell propagators are needed to build the complete ansatz. 



A generalized unitarity cut is a sum over products of amplitudes, 

iC Yl A (1) A {2) A {3) ■ ■ ■ A (m) , 



(2.17) 



states 



evaluated for kinematics that place all c cut lines on shell. Here we normalize the cuts 
to include the factor of i of each cut Feynman propagator. Each cut-line particle appears 
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(a) (b) (c) (d) 




FIG. 8: These 11 cuts, along with the two-particle cuts in fig. [9l suffice to determine any massless 
four-loop four-point amplitude. As displayed they determine any planar amplitude; to determine 
any non-planar contributions one must also include cuts with the legs of each tree sub-amplitude 
permuted arbitrarily. The cuts (h) and (k) have been evaluated in D dimensions, as they are 
included in the box cuts fig. fTUT c) and (e), respectively. 




(a) (b) 

FIG. 9: The two-particle cuts. These cuts have been evaluated in D dimensions. Together with 
the cuts of fig. [8] these form a spanning set of cuts for any massless four-loop four-point amplitude. 

twice in the summand — leaving one amplitude and entering another. Often the cuts are 
chosen so that all the amplitudes are tree amplitudes, although this is not necessary. Fig. [7] 
illustrates particularly useful cuts for determining the four-loop four-point amplitude, in 
which the maximal and near-maximal number of propagators are cut. (In order to avoid 
excessive clutter, we draw maximal and near-maximal cuts without the usual dashed lines 
to indicate the cuts.) Another set of useful cuts for confirming all contributions is depicted 
in figs. M and |9j 

For massless theories, the on-shell momentum condition for the particle associated with 
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each cut line is If = 0. We sum over all possible states of the cut-line particle. For M = 
4 sYM theory the sum is over the full supermultiplet, including gluons, scalars and fermions. 
For pure- Yang-Mills theory in four dimensions the sum is only over positive and negative 
helicity gluons. In the latter case, there can also be contributions, such as loops on massless 
external legs, that are not detectable in the cuts of figs. E] and While these contributions 
vanish after loop integration in dimensional regularization, they can nevertheless alter the 
UV behavior of the amplitude because their vanishing is the consequence of a cancellation 
between UV and IR divergences. In the case of M = 4 sYM theory, for representations of 



the amplitude that obey manifest power-counting |16l . |30| and do not contain three-point 
sub-amplitudes with one external and two internal legs (which is the case in this paper and 
in previous work), such massless-external-loop contributions do not appear, by virtue of 
unitarity and the vanishing of supersymmetric three-point loop amplitudes. 

The cut- construction of a multi-loop amplitude formally begins with a generic ansatz for 
the amplitude in terms of multi-loop Feynman integrals, in which the numerator polynomi- 
als of each integrand Ni(£j,kj) contain arbitrary coefficients. 1 These coefficients are then 
systematically determined by comparing the generalized cuts of the ansatz to the cuts of 
the amplitude. 

For a multi-loop expression to be correct for a given theory it must satisfy all possible 
generalized unitarity cuts. However, many cuts are simply special cases of other cuts. We 
define a spanning set of cuts as any set whose verification is sufficient to ensure that all other 
cuts are satisfied, and thus the correctness of a multi-loop expression for the amplitude. 
Below we shall describe such sets. 

In the process of cut-verification, an ansatz for the desired multi-loop amplitude — de- 
scribed in terms of Feynman integrals over loop-momenta — is evaluated using on-shell 
momenta on the cut lines at the integrand level and compared to eq. (12.171) . using the same 
cut kinematics. This procedure requires first identifying the subset of parent graphs in the 
ansatz that are nonzero on the specified cut, and then lining up momentum labels between 
each such parent graph and the generalized cut (12.171) . The latter step can be performed by 
decomposing each of the constituent amplitudes in eq. (12.171) into their own parent graphs. 



1 The number of independent Feynman integrals present in the four-point amplitude is typically far smaller 
than the number of Feynman diagrams. 
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These graphs come with a labeling, which can be used to provide a suitable labeling for each 
parent graph in the ansatz. (One may have to perform some initial relabeling in order to 
avoid duplicate labels coming from different constituent amplitudes.) Finally, one permutes 
the labels in the original representation of the ansatz so that they match this labeling. In 
appendix [A] we illustrate this procedure in the evaluation of a nontrivial cut at four loops. 

As mentioned in the discussion of the parent-graph decomposition, the numerators of par- 
ent graphs are not uniquely determined by cut constraints, because of an inherent ambiguity 
in the assignment of contact terms. This ambiguity corresponds to the ability to add terms 
to the integrands associated with different parent graphs, such that the complete amplitude 
remains unchanged. Essentially, this ambiguity is nothing more than the freedom to add 
zero to the amplitude in a nontrivial way. This freedom allows various representations of the 
amplitude, which can expose different properties. A particularly useful property is for every 
term in every parent integral to be no more UV divergent than the complete amplitude; that 
is, all UV cancellations are exposed. It is also rather desirable for the numerator factors to 
respect the symmetries of the parent graphs. Indeed, we will find a representation of the 
four-loop four-point amplitude with both these properties. Imposing these properties at the 
beginning helps to limit the number of possible terms in the ansatz. 

In practice, as will be discussed in section IIII D\ it is possible to construct a multi-loop 
amplitude iteratively - - building and refining an ansatz by requiring it to be consistent 
with a sequence of generalized cuts, beginning with the maximal cuts. Note that, a priori, 
one can make any simplifying assumptions that restrict the ansatz; these assumptions are 
justified a posteriori by verifying the ansatz against a spanning set of cuts. Figs. Eland [9] 
show generalized cuts used in the verification of the four-loop four-point amplitude. In these 
figures all exposed internal lines are cut. 

Color-stripped cuts involve products of the color-ordered tree-level partial amplitudes 
that appear in the trace-based color decomposition for gauge group SU(N C ). On the other 
hand, color-dressed cuts are products of full tree amplitudes, including their color factors, 
which are products of structure constants f abc and are valid for any gauge group. One can 
always verify a color-dressed amplitude for any gauge group G by considering a spanning 
set of color-dressed cuts. In practice, we used color-stripped cuts in our construction of the 
four-loop amplitude; in the final step we color-dressed our integrands with color factors Ci 
which are built from structure constants f abc and are valid for any gauge group. 
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Although in principle all cuts can be evaluated analytically, some cuts in fig. [8] are rather 
nontrivial. It is therefore often more practical to evaluate the cuts numerically to high 
precision for a number of points in the phase space satisfying the cut conditions. Because 
the check is performed at the level of the integrand, and does not require any numerical 
integration, it can be performed at arbitrarily high precision. 

In the rest of this section we will review technology for summing over the M = 4 multiplets 
of states crossing each cut, which is needed to evaluate the unitarity cuts. In section [TTTI we 
will present two classes of unitarity cuts — the two-particle cut, and the box-cut — which 
have been worked out in all generality in terms of lower-loop parent-dressings. These cuts 
are particularly handy because they can be applied with very little calculation. 



1. Su^er symmetric sums over states in four dimensions 

We can simplify the evaluation of many cuts by restricting the internal loop momenta 
(as well as the external momenta) to four dimensions. Then we can enumerate the internal 
states according to their four- dimensional helicity, and apply powerful supersymmetry Ward 
identities jso| or on-shell superspace formalisms (both of which are valid only in four dimen- 
sions), in order to simplify the sum over intermediate states. In M = 4 sYM, this sum runs 
over the M = 4 super- multiplet, so we refer to it as a "super-sum" . For simple cuts, the sum 



over supersymmetric states in eq. (12.1 7p is easy to evaluate component by component |16| . 
by making use of supersymmetry Ward identities that relate the different tree amplitudes, 
and hence relate the different terms in the state sum. 
As described in some detail in ref. 



64j, for maximal or near-maximal cuts, it turns out 
that one can avoid nontrivial sums over particles by making use of solutions to the cut 
conditions that force all, or nearly all, particles propagating in the loops to be gluons with 
a single helicity configuration. Such restrictions can be arranged for cuts that contain suf- 
ficiently many three-point tree amplitudes. These solutions were sufficient for constructing 
the complete four-point amplitude in M = 4 sYM in this paper. 2 Remarkably, this allows 
us to build a complete ansatz for the four-loop amplitude avoiding all nontrivial supersym- 
metric sums over particles crossing the cuts. Even so, we must satisfy all solutions to all cut 



2 This property is special to maximally supersymmetric amplitudes and will not hold in other theories such 
as QCD. 



20 



conditions, including those that impose no restriction on the particle content. As discussed 
earlier, we verify the correctness of our construction on a spanning set of cuts. In such cuts 
we must sum over all allowed configurations of particles crossing the cuts. A good means 
for summing over the states is therefore needed, especially given the nontrivial bookkeeping 
of states required at four loops. 

In supersymmetric theories, superspace provides an efficient way to track contributions 
from different states in the same super-multiplet. However, we prefer a superspace which 
works well with on-shell methods. Such a superspace is based on Nair's construction, which 



encodes the MHV tree amplitudes of M = 4 sYM 70] . In recent years, this superspace has 
been generalized to any four- dimensional tree amplitude and also to loop level |42|, l7ll-l76|. 



A solution to the problem of evaluating super-sums in generic multi-loop unitarity cuts was 
given based on an MHV- vertex generating-function approach. In ref. [iJ], this solution 
was recast into two complementary approaches for efficiently evaluating multi-loop unitarity 
cuts. In the first approach, the problem is recast into the calculation of the determinant of 
the matrix associated with a certain system of linear equations. In the second approach, 
used in this paper, the contributions of individual states are tracked via SU (4) "i?-symmetry 
index diagrams". 

To systematically step through the many cuts we used to verify our construction of the 
four-loop amplitude, it is helpful to have an efficient and easily programmable algorithm 
for evaluating any cut, with essentially no calculation. The i?-symmetry index diagram 
method [z^j] is based on the observation that, after applying the MHV- vertex expansion 



for tree amplitudes 8l], the cuts of M = 4 sYM amplitudes are simply related to those 
of (non-supersymmetric) pure Yang-Mills theory. By carrying out the super-sums in the 
MHV-vertex expansion, each term contains a numerator of the form, 

(S 1 + S 2 + --- + S m ) 4 , (2.18) 

where the S^s are spinor-product monomials, such as (i\ i 2 ) [13 14] ■ ■ ■ (ij-i ij)- Upon expan- 
sion of eq. (I2.18p . each quartic expression SiSjSkSi corresponds to a single assignment of 
helicities to particles crossing the cuts. Remarkably, eq. ( 12.181) can be inferred instead from 
the much simpler state sum for pure Yang-Mills theory, for which the analogous numerator 
is 

Sl + S* + .-- + Si. (2.19) 
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FIG. 10: The three box cuts. These cuts have been evaluated in D dimensions. 



One introduces anticommuting parameters, which transform under the SU (4) R symmetry 
of M = 4 sYM, and track the relative signs between Si and Sj in eq. (I2.18p . With the aid 
of these parameters, the result (12. 18)) for the M = 4 sYM cut can be read off from the pure- 
Yang-Mills cut (12.191) . A detailed description of the algorithm, as well as the i?-symmetry 
index diagrams, may be found in ref. 74 1. 



Using this algorithm we have evaluated cuts containing only MHV and MHV vertices, 
as well as the spanning set of all 13 cuts in figs. [8j and |9] and their permutations. These 
evaluations confirm that the ansatz we constructed for the amplitude, using graphical rules 
and information provided by the maximal and near-maximal cuts, captures all contributions 
that are nonzero when loop momenta are restricted to four dimensions. 



Elvang, Freedman and Kiermaier 



73] have used the MHV-vertex expansion to provide 



very compact expressions for the super-sums for cuts (a) and (j) in fig. [BJ which contain 



next-to-MHV (and next-to-MHV) amplitudes. We also compared the cuts of our ansatz to 
their results, and found agreement. 3 Cut (a) is particularly powerful because it checks most 
of the terms in all 50 of the parent graphs in eq. (I2.16p . providing an important independent 
check. 



2. D- dimensional cuts 

The method outlined above efficiently solves the problem of evaluating unitarity cuts in 
four dimensions in M = 4 sYM. However, because we are interested in computing the am- 
plitudes in (D > 4)-dimensions, this is not sufficient; the loop momenta are Z)-dimensional. 
Even if we are interested in amplitudes in four dimensions, we need to use a (supersymmetry- 



preserving 82] ) form of dimensional regularization to regulate infrared singularities. If the 



We thank H. Elvang, D. Freedman and M. Kiermaier for assistance in implementing their expressions. 
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cuts are evaluated in four dimensions, as described in the previous subsection, then terms 
that vanish when the loop momenta are restricted to four dimensions, but are non-vanishing 
in D dimensions, could be missed. Indeed, amplitudes in theories with fewer supersymme- 



tries 



S3 



841 ]. and M = 4 sYM amplitudes with more than four external legs 23J, do contain 



such terms. Unfortunately, generic D-dimensional cuts are significantly more complicated 
than their four- dimensional counterparts. For continuous values of D, there is no helicity 
formalism, nor is there a useful on-shell superspace. Some of the additional complexity can 
avoided by working in D dimensions, but with the M = 4 states organized according to the 
M = 1 super- Yang- Mills theory in ten dimensions 13J. Nevertheless, evaluating general cuts 
in this way at four loops is still difficult. 

It is important to note that we do not need the full power of D-dimensional cuts to 
construct the amplitude, but only to identify any potential terms dropped by the four- 
dimensional cuts. Such terms are rare or even nonexistent for M — 4 sYM at low orders 

oops, and for the 



15 



16 



29 



85| 



and low multiplicity. Specifically, for four-point amplitudes through three 
planar contributions through four loops, explicit computation has revealed [13 
that the D-dimensional versions of the amplitudes are obtained simply by replacing the four 
dimensional loop integration measure with the D-dimensional one, 

d 4 p f d D p 



(27T) 



(2tt 



(2.20) 



and reinterpreting all Lorentz products of momenta as D-dimensional ones. Based on this 
evidence, we have every reason to believe that eq. (I2.20p holds as well for the non-planar 
contributions at four loops. Although we have not checked all cuts in D dimensions in this 
paper, we have performed a set of strong consistency checks to make it extremely unlikely 
that any D-dimensional contributions have been missed. Such checks include all two-particle 
cuts, and all generalized cuts that isolate a four-particle sub-amplitude (box cuts), as shown 



in fig. [TUJ As noted some time ago 15|, |l6[, for four-point amplitudes, the iterated two- 



particle cuts automatically give the same result in D dimensions as in four dimensions. In 
the next section we will explain why the box cuts have the same property. Another powerful 
check comes from the new diagrammatic numerator identities 
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which hold in any 



dimension. They allow us to obtain many non-planar terms directly from planar ones. At 



four loops, the latter are known to be valid in D dimensions 



13], at least for external 



gluon states. Because of these checks it is rather unlikely that any terms were dropped in 
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extending the four- dimensional loop-momentum integrand to D dimensions. Nevertheless, 
it would still be useful to evaluate a complete set of unitarity cuts in D dimensions. As a 



step in this direction, the unitarity cuts have been confirmed for six-dimensiona 



and cut momenta 



86], using the helicity formalism of Cheung and O'Connell 87] and the 



external 



on-shell superspace of Dennen, Huang and Siegel 



III. CONSTRUCTING A COMPACT ANSATZ 



In the process of constructing amplitudes, it is helpful to have a toolkit that allows one 
to write down large classes of terms with essentially no computation. Even heuristic rules 
motivated by observed structures, or tools that capture only a subset of terms, can be 
rather useful. Typically, such rules allow one to quickly fix the structurally simplest terms 
in the amplitude, allowing the remaining effort to be focused on the more intricate ones. 
This strategy is especially potent when combined with the method of maximal cuts 64], 
which (as discussed below) allows a relatively small set of contributions to be considered in 
isolation. 

For planar M = 4 sYM there are a set of powerful graphical tools. The oldest of these 
tools is the "rung insertion rule" 15], which generates certain higher- loop contributions 
from lower-loop ones. More recently, the observed dual conformal properties of planar 
Af = A sYM M 



to prefactors 
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3] amplitudes have led to a powerful method for determining them, up 
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64, 



that can be determined straightforwardly from cuts. Heuristic 



ru 

65 



es 



or determining the prefactors in the planar four-point case have been given as well 



64, 



69] . Unfortunately, it is not clear how to extend the notion of dual conformal invariance 



to non-planar contributions. We also remark that the planar terms in the four-loop four- 
point amplitude were determined previously (without using dual conformal invariance) 13 ] . 
The non-planar contributions are much more intricate. In this section we will discuss tools 
that are useful for identifying both planar and (more importantly) non-planar contributions. 

We begin by reviewing and extending some particularly useful cuts that can be expressed 
very simply, and in all generality, in terms of lower-loop expressions. We also discuss a tree- 
level identity that allows many multi-loop contributions to be constructed, up to potential 
contact terms. We will close this section by discussing the method of maximal cuts, which 
provides a systematic tool for constructing all terms in any amplitude, including contact 
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terms. 



A. Two-particle cuts 



1 



2 




FIG. 11: A two-particle cut that may be used to construct contributions to the (L + L' + l)-loop 
amplitude from those at L and L' loop orders. 

A two-particle cut of a multi-loop four-point amplitude has the form shown in fig. [TT1 — 
it divides the amplitude into two lower-loop four-point amplitudes. Four-point amplitudes 
in M = 4 sYM have an especially simple dependence on the external states. This fact makes 
it possible to immediately write down the numerator factors for parent graphs that have 
two-particle cuts, in terms of lower-loop numerator factors (up to potential contact term 
ambiguities). This method can be applied to non-planar parent graphs as well, making it 
especially powerful. 

This simplicity relies on the observation that all M = 4 sYM four-point amplitudes can 
be expressed in a common factorized form, 



where A4 represents the full color-dressed amplitude (as distinguished from the color- 
stripped A4). All of the state dependence in A4 is carried by the kinematic prefactor 
/C = s 12 S23A t / [ cc , also defined in eq. (12.61) . All of the color dependence is carried by the 
state-independent universal factor . We will use this factorization of color and state 
dependence to determine the terms in that are visible in two-particle cuts, iteratively 
in terms of lower-loop universal factors. This result will be valid in D dimensions, whenever 
the lower-loop universal factors are valid in D dimensions. 

For four- dimensional external momenta and states, eq. (13.11) follows from the Ward iden- 
tities for maximal supersymmetry. These identities relate all four-point amplitudes to each 



4 L) (1,2,3,4) = £ 2L+2 /C(l,2,3,4) f/ (L) (l,2,3 



4), 



(3.1) 
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other at any loop order 16|, |80|, |9Q|] • From explicit computations, we know that this equation 



holds in D dimensions 
three loops 



10 



15 



16 



'or all states at one and two loops and at least for gluon amplitudes at 



85] . Using the observation that, in theories with 16 supercharges, the 



number of states (2 8 ) in a massive representation of the supersymmetry algebra is the same 
as the number of states in a product of two short (massless) representations (2 4 x 2 4 ), Alday 
and Maldacena argued that the intermediate states in a 2 — y 2 scattering process form 
a single supermultiplet in any dimension. This argument suggests that eq. (13. ip also holds 
in any dimension, with the loop factor capturing the L-loop correction to this multi- 
particle intermediate state. We therefore assume that eq. (13.1 p is valid in D dimensions, for 
any two-particle cut of the four-loop four-point M = 4 amplitude. 

In order to treat the tree-level case, L = 0, on an equal footing with loop level, we note 
that the color-dressed tree amplitude in M = 4 sYM can be written as 



A^\l, 2, 3, 4) = g 2 Af(l, 2, 3, 4) x ( f a ^ b f ba ^ + f<**i* jW 4 £^ 

V s 13 



(3.2) 



where A± = A^ cc , and we have used the color- Jacobi identity to eliminate the color factor 
j-a 1 a 2 bj?ba 3 a4, j n f avor f the other two. We also used the fact that K, is crossing symmetric 
(see eq. (I2.10p ). which implies that all the orderings of the color-ordered tree amplitude 
A\ (1, 2, 3, 4) are related simply to each other, up to ratios of kinematic invariants. Dividing 
eq. (I3.2p by /C, we see that the state-independent color-dressed universal factor at tree level, 
U(°\ defined by eq. (13. ip . is given by, 

(fa^aib -Cba,2a,3 £aza\b ibaia^ \ 
J - J - + I J - . (3.3) 

s 12 s 23 s 12 s 13 / 

In general, the universal factor is a sum of L-loop integrals. The integrands en- 
tering U^ L > are rational functions of momentum invariants involving the loop and external 
momenta. Explicit formulae for the universal factors for L = 1,2,3, including planar and 
non-planar contributions, may be found by matching eq. (13.11) with the known amplitudes 
already presented in section UTS] (E2~TB and (ETT4D . 

Next we evaluate the generic two-particle color-dressed cut depicted in fig. [TTJ It cuts 
the (L + L' + l)-loop amplitude A^ +L 2, 3, 4) into the two four-point amplitudes 
A^\— h, 1, 2, l 2 ) and a!\ \— / 2 , 3, 4, / x ) of loop orders L and L', respectively. The cut has 
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the form, 



A% +L ' +1 \l, 2, 3, 4) = i 2 J2 A { 4 L \-k, 1, 2, l 2 ) A?'\-l 2 , 3, 4, h) 

2-cut z — ' 



(3.4) 



N=4 
states 



where the state sum is over the particles with momenta l± and l 2 . 

Using the factorization (13. ip and the state-independence of U^ L \ we can immediately 
rewrite the cut as follows: 

tf(i+£'+i)( lj2> 3,4)_ x £(1,2,3,4) = z 2 [/( L )(-/ 1 ,l,2,/ 2 )[/( L ')(-/ 2 ,3,4,/ 1 ) 

K(-h, 1, 2, / 2 ) £(-Z 2 , 3, 4, /x) . (3.5) 



2-cut 



X 



Af=4 
states 



Substituting in the definition of /C given above, we find: 

U^ L+L ' +1 \l,2,3,A) x Sl2S23 4 0) (l,2,3,4) = i 2 U^(~h, 1, 2, l 2 ) U^ L '\-l 2 , 3, 4, h) 



2-cut 



^12^^ Yl Af\-h,l,2,l 2 )Af\-l 2 ,3A,h). 

(3.6) 



JV=4 
states 



To evaluate this, we use the sewing relation between two four-point color-ordered N = 4 sYM 
trees 



15 



16] 



4°V^l,2,Z 2 )4 0) H2,3,4,/i) = -is 12 s 23 Af\l, 2,3,4) — 



Af=4 
states 



S2l 2 s 4:l 1 



(3.7) 



This sewing relation is valid in any dimension D and for any external states in the M = 4 
multiplet. A straightforward way to confirm eq. ( 13. 7ft is to work in D = 10 and evaluate the 
sum over states in components, using the fact that in D = 10 M = 4 sYM is equivalent to 
an M = 1 theory composed of a gluon and a gluino. By dimensional reduction the sewing 
relation ( 13. 71) then holds in any dimension D < 10. Recently, this equation has also been 
verified directly in six dimensions using an on-shell superspace 88]. 

Applying eq. (13. TP to eq. ( 13. 6ft . we find the key equation for building all contributions 
from two-particle cuts directly in terms of the U s: 



tfL+V+l)^ 2j 3j 4) = • s 2 2 ^(L) ^ 2) h) U(L') { _^ 3> 4j ^ 



2-cut 



(3.8) 



Equation (13. 8p is rather powerful. No complicated calculations remain in order to obtain 
all contributions visible in two-particle cuts; they are given simply by taking the product of 
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FIG. 12: Sample contributions to the full color-dressed two-particle cut for L = 2 and L' = 1. 
The diagrams on the right-hand side show some of the terms in LT( 4 )(l, 2, 3,4) that are constructed 
from these cuts, using fig. [T3j The explicit color factors, as well as factors of i, have been omitted. 
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FIG. 13: The lower-loop integral functions entering the cuts on the left-hand side of fig. [12j 
This figure displays in detail how the prefactor of the planar double-box integral appears in 
U&X-h, 1,2, l 2 ) as either S12 or (fci — h) 2 , depending on the permutation. 

lower-loop results. The color-dressed JJ^ L+L ' +l ^ is given immediately as a sum over products 
of individual integrals residing inside the and factors, up to terms that vanish 

because of the on-shell conditions, l\ = l\ = 0. (As a straightforward exercise, one can verify 
that the one-loop universal factor f/W — which can be extracted from eq. (I2.5P — satisfies 
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this equation, using the tree-level universal factor given in eq. (13. 3p .) 

Figs. [12] and [13] illustrate diagrammatically some of the terms generated by eq. fl3.8[) 
for the case L = 2 and V = 1. For simplicity, we draw only the planar contributions of 
^ (— /i, 1, 2, Z 2 ), encoding the f abc visually in the diagrams, and we omit all factors of i. The 
denominator factors in and correspond to propagators that are visible on the left- 
and right-hand sides of fig. [13] respectively. Therefore they are accounted for graphically in 
| 2 - cu t simply by connecting the l\ and Z 2 legs of the corresponding diagrams. Similarly, the 
numerator factor for each parent graph on the right-hand side of fig. [12] is given by forming 
the product of the numerator factors for the two sewn subdiagrams in fig. [T3] (taking into 
account the proper permutation of legs), and then multiplying by two powers of Si 2 . 

This diagrammatic interpretation of the two-particle cuts provides a rather simple tool for 
generating many higher-loop contributions from known lower- loop ones. It is the mechanism 



behind the rung rule [15|, [16| . For the planar contributions at four loops, the two-particle 
cuts have either L = 2, V — 1, as in fig. CE2] or else L — 3, L' = 0. Together, they capture 
diagrams (a)-(e) in fig. [5] but not diagram (f). (Diagram (f) can still be guessed from the 
rung rule, or constructed using a box cut, as described in the next subsection.) These cuts 
also do not guarantee the absence of contact terms that have no two-particle cuts, such as 
diagrams (f 2 ) and (d 2 ) in fig. [6] For the full four-loop amplitude described in section HVl 33 of 
the 50 parent graphs contain two-particle cuts (graphs 1 through 27, and graphs 40 through 
45). The two-particle cuts capture the majority of the terms contributing to these graphs. 
Because the two-particle cut sewing algebra is valid in D dimensions, all contributions 
obtained by iterating two-particle cuts are automatically valid in D dimensions. Surprisingly, 
the two-particle cuts capture the majority of terms in the 33 parent graphs containing them. 
The fact that so many potential contact terms are absent hints at further structures to be 
uncovered. 

We note that in M = 8 supergravity, the two-particle cuts have an equally simple struc- 



ture [la], which can be exploited analogously. 



B. Box cuts 

The simple structure of the four-point amplitude in M = 4 sYM can also be applied to 
generalized) four-particle cuts that isolate a four-point sub-amplitude. A simple version 
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FIG. 14: The box-substitution rule [641 ] for generating a higher-loop contribution by inserting a 
one-loop four-point box subintegral into a four-point vertex. In this example, we substitute a box 
into the central four-point vertex in the four-loop "window" diagram. The result is a five-loop 
integral that cannot be obtained from two-particle cuts. (Note that an overall normalization factor 



of S12S23 has been absorbed into fC, relative to ref. 



6J. 



of this generalization appeared already 64J as a "box-substitution rule". It allowed the 
construction of L-loop contributions with a box subgraph, starting from (L — l)-loop contri- 
butions with a contact interaction, as illustrated in fig. HU Related rules were discussed in 
conjunction with leading singularities 69]. Here we promote the box-substitution rule into 
a more general cut for M = 4 sYM amplitudes in D dimensions, which we call the "box 
cut" . 

Consider the generalized cut of an L-loop n-point amplitude, 



box cut 



(i) 



■A (n ...A, > 



(3.9) 



JV=4 
states 



that is composed of a generic set of color-dressed amplitude factors, except for the i such 
factor, which we take to be a color-dressed L'-loop four-point sub-amplitude, A\ ,L . (There 
may be additional cut conditions imposed on this sub-amplitude; its internal kinematics are 
irrelevant for the subsequent discussion.) Example of such box cuts are given in fig. [T5l 

The L'-loop four-point sub-amplitude of M = 4 sYM is special because of the factorization 
property (13. ip . Labeling the cut legs by l\, l 2 , h, U, we have, 



A^ih, k, k, k) = AfUh, h, h, k) (h + k) 2 (k + h) 2 uM(h, l 2 , h, k) 



%(0 



(3.10) 



where, as in the previous section, we use A4 to represent the color-dressed amplitude, and 
only the color-ordered tree amplitude factor A^L, depends on the states crossing the cuts. 
Therefore we can pull the factor U^') out of the sum over states in eq. (13. 9p . leading to a 
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FIG. 15: Two examples of multi-particle multi-loop "box cuts". They reduce to lower-loop cuts by 
replacing the darker (red) four-point sub-amplitudes by four-point color-ordered trees, multiplied 
by known numerator and denominator factors. This property allows these cuts to be computed 
easily in D dimensions, once the amplitudes with fewer loops are known. White holes represent 
loops and the darker sub-amplitudes mark four-point amplitudes amenable to reduction. 

simpler expression in the summand, 

A^ = (h + l 2 ) 2 (l 2 + h) 2 U^ L '\h, l 2 , l 3 , h) V A (1) ■ ■ ■ Af ■ ■ ■ A {m) . (3.11) 

box cut , ,v 1 



states 



The state-sum is identical to a lower-loop cut, that of the (L — L')-loop amplitude, but 
utilizing the color-ordered contribution to the i th tree. This fact immediately gives a simple 
relation between the L-loop box cut and contributions to the reduced (L — L')-loop cut 
under the same cut conditions. 

We can formally write down an equation relating the cut of an L-loop amplitude to a cut 
of a lower-loop one as, 

A^ = (h + k) 2 (k + h) 2 U^ L '\h, l 2 , l 3 , h) ■ (3.12) 

box cut cut 

We introduced the reduced cut A notation to emphasize that the state-sum in eq. (13. lip is 
exactly a (L — V) loop unitarity cut which is color-dressed with f abc everywhere, except for 
the four-point color-ordered tree amplitude whose associated color factors are accounted for 
in the L'-loop universal factor U^ L '\ 

Given a generalized cut that isolates an L'-loop four-point sub-amplitude with legs 
h,h,h,h, we can re-express the box cut as a recipe that can be applied easily to indi- 
vidual diagrammatic (integral) contributions: 



Split up the cut into three parts as in eq. (13.121) : The reduced cut, A^t L ''| cut ) the 
kinematic factor (Zi + hYih + h) 2 , and the loop integrals £/( L ') of the four-point sub- 
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amplitude. (This latter part generalizes to 11 loops the one-loop box integral of the 
box substitution rule.) 

• Express the reduced cut of the known lower-loop amplitude in a diagrammatic form 
that corresponds to a covariant integral representation. 

• The diagrams of the reduced cut may contain spurious propagators in the (Zi + I2) 2 
or {I2 + I3) 2 channels, which upon multiplication cancel against the (h + h^ih + h) 2 
prefactor in eq. (I3.12p . The result is always a diagram with an internal four-point 
contact vertex. 

• To recover the integrals of the original box cut, insert the four-point integrals of 
(e.g., the box integral for V — 1) into the obtained four-point contact vertex of each 
diagram. 

Fig. [16] shows how the box cut can be used to determine the numerator polynomial for 
fig- [5]^f ) , using the three-loop information in fig. IHe). Although this example is planar (and 
is presented in a color-ordered way in the figure), it is just as simple to use the box cut for 
non-planar contributions. For example, inserting a box into the four-point vertex in fig. [TS] 
in a non-planar fashion generates contributions to parent graph 29 in the full four-loop 
amplitude. 

The box cut is an extremely efficient way to obtain contributions to parent graphs that 
contain a lower-loop four-point subgraph. As mentioned earlier, the V = 1 box cut is 
closely related to the box-substitution rule. The box cut also generates contributions that 



are consistent with the rung rule 



15|- 



Box cuts capture a majority of those terms in the complete four-loop four-point amplitude 
in section HVl that are not determined by two-particle cuts. Of the 17 parent graphs that do 
not have two-particle cuts, 13 of them have box cuts. The only four that have neither two- 
particle cuts nor box cuts are graphs 39, 48, 49 and 50. In fact, most of the parent integrals 
have multiple box cuts, allowing us to constrain their numerators under complementary cut 
conditions, and to fix many of the contact terms. 

From the above covariant derivation it follows that box cuts are valid in any dimension, 
if both the reduced cut A^n~ L ^ I , and the four-point universal factor are known in D 

I cut 1 

dimensions. As a practical matter, the universal factors entering the lower-loop amplitudes 
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FIG. 16: Application of the box cut to determine the numerator N for the four-loop parent graph 
in fig. E](f). The (red) dashed cut conditions around the upper left box in this diagram allow us to 
replace it by the product of a reduced cut diagram, some kinematical factors and a box integral. The 
reduced cut diagram is then expanded into two three-loop "tennis-court" diagrams, corresponding 
to the two allowed channels of the marked four-point vertex. The relevant kinematical pieces of 
the tennis-court diagrams, i.e. the numerators and the spurious propagators, are extracted from 
the known three- loop contribution in fig. [U[e). Assembling all the kinematical factors gives the 
result for N, which is free of spurious propagators. The result is consistent with fig. E^f); here the 
numerator is symmetrized with respect to the (1 <->■ 3) symmetry of this parent graph. 

should already be known in D dimensions, prior to attempting the higher-loop calculation in 
D dimensions. In the case relevant to this paper, L = 4 and n = 4, all we need as input are 



the L' = 1, 2, 3 four-point amplitudes, which are indeed known in D dimensions [10|, [15|, |16 |. 

The effectiveness of the box cut suggests that one should investigate analogous "pentagon 
cuts" , etc., which isolate sub- amplitudes with five or more legs. Both the color and kinematic 
structure of the five- and higher-point loop amplitudes is, however, more intricate, and there 
is no simple factorization property similar t o eq . ( T3. 10[) . (See for example, the five- and six- 



point loop amplitudes described in refs. [23l. 184 l89|.) 



The box cut also easily generalizes to M = 8 supergravity, because its four-point ampli- 
tude has a factorized form similar to eq. (13.11) , which is related to the existence of analogous 
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supersymmetric Ward identities 



80| 



C. Color-kinematic duality 



In the early 1980s, radiation zeroes appearing in certain gauge-theory cross sections were 
traced back to a curious identity obeyed by tree-level four-point amplitudes 
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9j. This 



curiousity turns out to be the simplest of a set of relations arising from a general tree-level 
duality between color factors and kinematic numerators 36|. If one assumes that the duality 



holds for an arbitrary number of external states, one can derive 
color-ordered tree amplitudes, which have since been proven 
string theory amplitudes have also been proven recently 
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new relations among 



93(. Similar relations among 



951 ] . In the low-energy limit, the 



string-theory re 
Kuijf relations 



ations become identical to two types of field-theory relations: the Kleiss- 



96| (which follow from color considerations alone 97|) and the amplitude 



relations which follow from the color-kinematic duality 

In this subsection, we discuss how the four-point tree-level color-kinematic identity may 
be combined with generalized unitarity at the loop level 36| , particularly to the construction 
of the four- loop M = 4 sYM amplitude. In short, the identity relates sets of three parent 
graphs that only differ in how a four-point cubic tree graph is glued into the rest of the graph. 
Evidence that the color-kinematic duality also holds directly at the loop level, without the 
need to impose on-shell conditions, was presented recently for the three-loop four-point 



amplitude 
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N, 




FIG. 17: Graphs for the four-point tree amplitude. Contact terms are absorbed into the diagrams 
as inverse propagators. Each diagram is associated with a color factor obtained by dressing the 
vertices with an f abc , as in section Hi Al 



Consider the color-dressed four-point tree amplitude. Just as for the loop amplitudes dis- 
cussed in section III Al it can be written as a sum of color factors Cj multiplied by kinematic 



34 



factors. The kinematic factors can be further divided into denominators, which are prop- 
agators associated with (tree- level) parent graphs, and numerators N. As in section III A\ 
contact terms can be absorbed into the Ni, so that we require only the three cubic graphs 
shown in fig. [T71 In this representation the amplitude is, 

^ r = 9 W + ^ + ^A (3.13) 
\ s t u J 

where s — (k± + k<i) 2 , t = (/^ + k 3 ) 2 and u — (fci + k 3 ) 2 correspond to the three channels, 
and 

Q — Jai_a,2b jtbasai _ Ja2as,b Jba^ai q _ ja^aib jba-sa\ ^-^) 

are color factors corresponding to the three graphs in fig. [T71 The color factors of the graphs 
satisfy the Jacobi identity, 

C u = C s -C t . (3.15) 

The Ni in eq. ( 13 . 1 3[) contain momentum invariants, polarization vectors, spinors and 
superspace Grassmann parameters. The only real restriction on them is that eq. ( 13 . 1 3[) 
gives the correct color-dressed tree amplitude. Hence there is a tremendous amount of 
freedom in the definition of the numerator factors. (Non-local Ni could even be allowed.) 
This freedom is just the tree- level analog of the inherent ambiguity in the multi-loop parent- 
graph decomposition mentioned in section III Bl We refer to the invariance of eq. ( 13 . 1 3[) 
under this freedom as a "generalized gauge invariance." For every such generalized gauge 
choice for the four-point J\f — 4 sYM amplitude, the numerator factors must satisfy the 



identity [36] 



N u = N s - N t , (3.16) 

in concordance with the color Jacobi identity ( 13 . 1 5 j) . We emphasize that the identity ( 13 . 1 6[) 
is only between the numerator factors; it does not involve the propagators associated with 
the 5, t, and u channel graphs. It is fairly straightforward to check that these identities hold 



in D dimensions by direct computation 



9l| . Although it is not relevant to this section, it 



should be noted that for higher-point tree amplitudes the color-kinematic duality is only 



36, 



63] 



manifest for certain special generalized gauge choices 

In conjunction with the unitarity method, the tree- level four-point numerator iden- 
tity (13 . 1 6[) becomes quite powerful. In every multi-loop parent graph that contains four 
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FIG. 18: Two nontrivial relations between numerators of planar and non-planar contributions at 
four loops, which follow from the tree-level numerator relations in fig. [T71 For the three graphs in 
each relation, the configuration of lines outside the region marked by a dashed circle is identical; 
the only difference is for the lines fully inside this region. The relation on the second line involves 
the same graph numerator N28 with two different labelings of momenta. 



on-shell propagators arrayed around a four-point tree sub-amplitude, it relates the numer 



ator factor to those of two other parent graphs satisfying those conditions [36|. The three 
multi-loop parent graphs correspond to gluing in the four-point cubic tree graph in its s, 
t, or u channel configuration. For every line of each cubic graph this identity will always 
relate the numerators of three graphs. However, the relations do not have to be manifest in 
a given amplitude representation, because of the freedom to move contact terms 4 associated 
with other propagators between different graphs. 

These relations allow one to take kinematic numerator information, obtained using dual- 
conformal symmetry (for planar graphs), two-particle cuts and box cuts, and export that 
information to other parent graphs or contributions for which such methods are not appli- 
cable. 

To see how this works, consider the two four-loop examples illustrated in fig. [18j In each 
case, the numerators on the right-hand side are planar, and are relatively simple to obtain 



4 One can automatically disregard such contact terms by considering near-maximal cuts where only the 
central propagator in the four-point tree graph is off shell. 
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using two-particle cuts (Nq and Nis) or box cuts (N28, at least up to contact terms). We 
then use the four-point tree color-kinematic duality to obtain the bulk of the non-planar 
graph numerators on the left-hand side of each equation. The numerators of the four-point 
tree amplitudes entering the cut satisfy the relation (I3.16P (see also fig. IT?]) . The remaining 
contributions from the outside of the dashed circle are identical in all three contributions. 
Hence we obtain a numerator identity for the loop integrands, 

^L = (^-iV c )| cut , (3.17) 

where the N x are the kinematic numerators of the integrals corresponding to the graphs in 
fig. [THJ This relation is valid in D dimensions. However, it holds only for the numerator 
terms that are nonvanishing under the imposed on-shell conditions, in which the four legs 
crossing the dashed circle are put on shell. Also, it should be realized that some contact 
terms may be distributed for convenience into other graphs. That is, there are a large 
number of coupled equations obtained from these constraints, and it is not necessary to 
satisfy each one simultaneously for all contact terms. 

To illustrate these ideas in more detail, consider the identity on the first row of fig. HE1 It 
involves the numerators N25, Nq and Nig which are presented in the next section, in figs. [20], 
[2T1 and [22] and also in appendix O We relabel the lines to match the labels in fig. [TBI 
obtaining, 

N 2§ = --512(545538 - S12S37) + /9S12S45 + /^8 S 12 , 



Arfig-LLOl _ 2 _ I _ _2 . 

1\ 6 - S 12 S 45 | { 3^ 4)5 ^ 7} - S 12 S 37 

rfig-P 



^18 — s 12S35S 46 ^ 3 ^ 4 6 _^ g j — S 12 S4 5 S38 . (3.18) 

a [TBI 

In addition, numerator iV 25 — picks up a sign relative to N25 in fig. [22] That is because 
the deformation of graph 25 in fig. [22] into the graph in fig. [18] requires an odd number of 
three- vertex reorderings (three). Each reordering results in a minus sign (from the structure 
constants) for the color factor C25, and a corresponding minus sign for N25. 

Using eq. f|3.18p . it is easy to see that the numerator relation (almost) holds on the cut: 



^•^L = (N^- . (3.19) 



The on-shell conditions on the legs crossing the dashed circle include /g = 0, so the term 
^■512^45 in eq. (13 . 1 8[) should be set to zero. What about the term /g/|si 2 ? It is not zero on the 
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cut, so it should be accounted for. The alert reader will notice that canceling propagators 
6 and 8 in graph 25 in fig. [22] gives a graph that is topologically identical to that obtained 
by canceling propagators 5 and 8 (or 6 and 7) in graph 28. Also, terms containing /|/| and 
Qlj are present in N 2 s- These features allow the l^s^ contact term in N25 to be moved 
elsewhere to be consistent with the identity. However, the presence of overlapping identities 
can complicate their application, when all contact terms are retained. 

The second relation in fig. [18] works similarly. The same graph 28 appears twice on 
the right-hand side, with two different labelings. It is worth noting that for our choice of 
numerators N 2 s and N 50 , as given in figs. [22] and [25] and appendix O, this particular relation 
holds even including all contact terms. 



It has been conjectured recently [68( that a representation exists for all multi-loop ampli- 
tudes in which all color-kinematic duality relations are manifest for all graphs, and with no 
internal on-shell conditions imposed. This conjecture has been confirmed for the three-loop 



four-point amplitude of M = 4 sYM, as well as for certain lower- loop cases 68] , but it 
remains to be tested more generally. Strong evidence in favor of the conjecture would be 
provided if the four-loop amplitude presented here can be rearranged into such a duality- 
satisfying form. We leave this exercise to future work. 

At four loops, the three rules just presented can be used to generate all non-contact- 
term contributions to the M = 4 sYM amplitude, as well as many of the contact-term 
contributions. To ensure that all contact terms are captured correctly, we turn to the 
method of maximal cuts. 



D. Method of Maximal Cuts 



The method of maximal cuts 29|, [64| offers a particularly efficient means for determining 
the numerator polynomials for each parent integral. In this method we start from generalized 
cuts with the maximum number of cut propagators (maximal cuts) and match these cuts 
against an initial ansatz. If an ansatz has been constructed that covers all non-contact-term 
contributions (for example, by using the three rules just presented), then this step is merely 
one of cut- verification. Next we systematically reduce the number of cut propagators (by one 
at each step) and match these (near-maximal) cuts — capturing in the process all potential 
contact contributions. 
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It is important that massless on-shell three-point amplitudes are non- vanishing and non 



singular 



63|, for appropriate choices of complex cut loop momenta [l^, |98|]. The maximal 



cuts of four-point amplitudes involve products of only three-point tree amplitudes, and are 
the simplest cuts to evaluate. Near-maximal cuts, in which one or two of the maximal-cut 
propagators have been allowed to go off shell, are the next simplest to evaluate, and so on. 

The advantage of the maximal-cut method is that it allows one to focus on a small 
number of terms at a time, namely those that become nonvanishing when a particular 
propagator is allowed to go off shell. This feature reduces the computational complexity at 
each stage, allowing us to efficiently find compact representations of amplitudes with the 
desired properties. We note that the "leading-singularity" technique, which is applicable to 
maximally supersymmetric amplitudes, is also based on cutting a maximal or near-maximal 
number of propagators 



69, 



991 ] . but in addition it makes use of further conditions from 



hidden singularities that are special to four dimensions. 

In practice the method of maximal cuts allows the sequential improvement of an ansatz for 
the numerator factor of each parent graph. Every new cut identifies the presence of missing 
pieces, which were left undetermined by the previous cuts, and which can be assigned to one 
of the parent graphs contributing to the cut. Because these pieces vanish on the previous set 
of cuts, they will contain an inverse propagator factor associated with the last propagator 
to be allowed to go off shell. Once new cuts cease to reveal any more missing pieces, the 
ansatz is generally complete and is ready for systematic cut- verification. 

Although the maximal-cut method can be applied to D- dimensional cuts, in order to 
simplify their evaluation we restrict many of the cuts to have four- dimensional momenta for 
both internal and external lines. As we often evaluate these cuts numerically, it is useful 
to build an ansatz for any missing pieces, which consists of a Lorentz-covariant numerator 
polynomial containing unknown constant coefficients. We reduce the number of unknowns in 
the ansatz by assuming that no individual term in it violates the expected ultraviolet power- 
counting bound (eq. ( 15 .11) below) [16J, |30|. These assumptions are, of course, validated by 
comparing against a spanning set of cuts after the amplitude has been constructed. 

At four loops, the bound (15. ip predicts that at most four powers of loop momenta (or 
at most two inverse propagator factors If) can appear in any numerator polynomial. This 
restriction allow us to focus our attention on the maximal and near-maximal cuts that 
have at least 11 cut conditions, If = 0, out of the maximal 13 (corresponding to the 13 
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FIG. 19: Four-loop examples of "MHV/MHV-amplitude cuts", which are composed entirely of 
four- and five-point tree amplitudes. 

propagators of the parent graphs). Examples of such cuts are shown in fig. [7J At the level 
of 11 cut conditions, there are always some quartic monomials of the form Ify^ that are 
non-vanishing. As one cycles through all cuts at this level, all such quartic terms will be 
detected, and their coefficients will be fixed. Similarly, one can show that these cuts will 
detect all quartic monomials of the form p 2 q 2 and more generally (p ■ p') (q ■ q') , where p,p',q 
and q' are linear combinations of the loop momenta and external momenta. We can continue 
the procedure of removing on-shell conditions, one by one, until we end up with a spanning 
set of cuts. However, in practice, it is much simpler to stop the construction phase as soon 
as we suspect that the ansatz is complete. 

The ansatz is then confirmed by checking that it matches the minimal spanning set of 
11 cuts in fig. [HI plus the two two-particle cuts in fig. [9j We refer to this set as a spanning 
set because the information it provides is equivalent to that contained in all possible cuts, 
and minimal because any further reduction could only involve tadpole-like contributions. 
To show that it is a spanning set, we show that it includes all the information in the 
ordinary two-, three-, four- and five-particle cuts. First of all, Fig. [8(a) is just the ordinary 
five-particle cut. The information from the ordinary two-particle cuts is given by fig. [9(a) 
and (b). Ordinary four-particle cuts consist of a tree-level six-point amplitude multiplied 
by a one-loop six-point amplitude. We can reproduce the information in these cuts by 
studying those generalized cuts in which we further cut the one-loop six-point amplitude in 
all inequivalent ways (omitting three-point trees). This procedure leads to fig. [8(b), (c), (d) 
and (e). Finally, ordinary three-particle cuts leave either the product of a tree-level five- 
point amplitude and a two-loop five-point amplitude (with further cuts leading to fig. E(f), 
(g), (h), (i), (j) and (k)), or the product of two one- loop five-point amplitudes (which does 
not lead to any new cut). In this classification, we can omit a cut if another cut already 
appears with a subset of the cut propagators. 
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Suppose one assumes the existence of a representation of the four-loop amplitude in which 
each term in the numerator polynomial for each parent graph has no more than two inverse 
propagators, consistent with the known J\f = 4 sYM power-counting fltl [3o| . In this case, 
one only needs to check near-maximal cuts with at most two canceled propagators. Because 
of this, the spanning set of cuts can be restricted to products of four- and five-point tree 
amplitudes, as illustrated in fig. [T9j (A six-point tree amplitude requires three propagators 



to be canceled from a maximal cut.) We refer to these as "MHV/MHV- amplitude cuts", 



because all tree amplitudes appearing in the cuts are either MHV or conjugate MHV am- 



plitudes. The MHV/MHV-amplitude cuts are useful because they are simpler to evaluate 
than the spanning set in figs. El and [HI the super-sums are particularly easy to evaluate jz^j]. 

We note that when the cuts are verified using color-stripped amplitudes, in order to 
capture all non-planar contributions we must include cuts where the legs of each tree entering 
the cuts are permuted in all possible inequivalent ways. 

IV. THE COMPLETE FOUR-LOOP AMPLITUDE 

We applied the construction methods outlined in the previous section to the four-loop 
four-point J\f = 4 sYM amplitude. The resulting amplitude is given by, 



5 4 

+ 4X10 + 2^11 + 4X12 + 9X13 + ^Xi4 + jXi5 + Xi6 + ^Xi7 + Xis + X19 



X20 + X21 + X22 + X23 + |X 2 4 + X25 + X 2 6 + |X 2 7 + 2^28 + X29 + |X 30 



+|x 31 + x 32 + x 33 + |x 34 + x 35 + x 36 + |x 37 + |x 38 + |X 3 g + |x 40 



+ |Xn + X42 + X43 + 2X44 + JX45 + |X46 + IX47 + |X48 + IX49 + IX50 , (4.1) 

where the prefactor /C, defined in eq. (I2.6p . encodes the full external-state dependence, 
and Xj = CJi are the color-dressed four-loop integrals. The Ci are color factors obtained 
by dressing the parent graphs with structure constants f abc , and are given explicitly in 
appendix[B] The /j(si2, S23) are D- dimensional loop integrals, defined in terms of numerator 
factors Ni in eq. f)2.4p . and corresponding to the 50 four-loop cubic parent graphs in figs. liZO} 
1251 The 50 th graph, which appears in eq. (14.1 ft and in fig. [251 is needed to match all cuts; 
however, it integrates to zero, and its associated color factor C50 also vanishes. Thus its 
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contribution to the integrated color-dressed amplitude is doubly vanishing. (Another reason 
we list this Af = 4 sYM contribution is because it gives a nonvanishing input into the 



construction of the corresponding Af = 8 supergravity amplitude [11].) 

As in the lower-loop amplitudes in section III A[ the sum runs over the 24 independent 
permutations of legs {1, 2, 3, 4}, denoted by S4, which act on both kinematic and color labels. 
The numerical coefficients in front of the integrals in eq. (14. ip are symmetry factors 1/S, 
where S is the number of elements in the discrete automorphism group of the corresponding 
unlabeled graph. As before, these factors compensate for overcounting. 

As mentioned in section HTT| the parent graphs containing two-particle cuts, namely graphs 
1-27 and 40-45, are the simplest to obtain; the bulk of the terms in their numerator poly- 
nomials are constructed using eq. (]3.8p . The remaining terms in these graphs, and all of the 
terms in the remaining 17 parent graphs, are obtained using box cuts (eq. (I3.12p ). plus the 
color-kinematic duality relation (I3.17p . as well as an evaluation of the near-maximal cuts. 
The amplitude's construction was followed by a confirmation of the complete set of cuts in 
figs. E and El 

V. ULTRAVIOLET PROPERTIES IN HIGHER DIMENSIONS 

In this section we examine the UV behavior of the four-point Af = 4 sYM amplitude 
in its critical dimension, i.e., the lowest dimension in which it diverges. This issue is of 
some importance because it provides a simpler venue than Af = 8 supergravity for studying 
UV cancellations in theories with maximal supersymmetry. For Af = 4 sYM, an analysis 
of supersymmetric cancellations in two-particle cuts (and more generally, in ordinary m- 
particle cuts that have only MHV amplitudes on either side of the cut), suggested the UV 



finiteness bound 16] 



D<D c = A+j (L>1). (5.1) 

One- loop amplitudes, corresponding to L = 1, are special as they are UV finite for D < 8, 
not D < 10. The bound (15.11) is somewhat stronger than earlier superspace power counting 



bounds 



lOOj, although all bounds agree that the theory is UV finite in D = 4. The bound 



(15. ip is consistent with a formulation of the theory having 3/4 of the superchar ges manifestly 



realized, and it has been confirmed 



using Af = 3 harmonic superspace [3l|. Explicit 



computations (including the ones discussed in this paper) demonstrate that this bound is 
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FIG. 20: Integrals (l)-(ll) appearing in the four-loop amplitude. The graphs encode denominator 
factors as the Lorentz square of the momenta flowing across every internal line, and color factors 
following the rules defined in fig. [TJ The momentum-dependent factor in front of each graph 
represents the numerator factor TV, that resides inside the integral ij, where (i) is the label below 
each graph. Numbers 1,2,3,4 label the external (outgoing) momenta. The internal legs carry 
momenta as signified by the arrows (here only leg 5 is labeled). The kinematic variables are 



defined as s 



(h + lj) and in the following figures s i = (h — lj) 2 and T%j = 2Z, • lj, where 



U is the momentum of leg i. A specific (clockwise) orientation of each cubic vertex (in the plane 
of the figure) is implied here. Due to the antisymmetry of the structure constants, any noncyclic 
reordering of a vertex should be accompanied by a sign flip of the numerator factor. 



saturated through at least four loops 



da 



,,151, 



16| . It is straightforward to verify from the 
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FIG. 21: Integrals (12)-(23) appearing in the four-loop amplitude. The notation follows that of 

fig. EDI 



planar amplitude in ref. 64J that the same result holds in the large- N c limit at five loops. 

An interesting question concerns the UV behavior of the different color terms. For an 
SU (N c ) gauge group it is convenient to expand the amplitude in the trace basis — i. e. traces 
of products of generators in the fundamental representation, as in eq. ( 12. lft . but including 
also terms that are subleading in the 1/N C expansion. For a four-particle amplitude there 
are only single- and double-trace terms. This decomposition follows quite closely that of the 
terms in the effective action that have four or fewer field strength tensors. The single-trace 
terms satisfy the finiteness bound (15. ip . As noted already in refs. 32j, double-trace terms at 
three and four loops exhibit additional cancellations which increase their critical dimension. 
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FIG. 22: Integrals (24)-(31) appearing in the four-loop amplitude. The notation follows that of 

fig. m 



These cancellations have been discussed in refs. [12J, |33|, which suggest double-trace terms 
should instead satisfy the finiteness bound, 



D < 4 



(L > 2) . 



(5.2) 



This formula is equivalent to the statement that an additional momentum invariant (s, t 
or u) can be extracted from the double-trace terms than the single-trace terms, at each 
loop order. As a statement about counterterms, it implies that the leading double-trace 
counterterms for L > 2 would have four covariant derivatives, for example Tr(D 4 F 2 ) Tr(F 2 ), 
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FIG. 23: Integrals (32)-(39) appearing in the four-loop amplitude. The notation follows that of 
fig- 1 



in contrast to the two derivatives characteristic of the single-trace counterterms Ti(V 2 F 4 ). 
In a superspace description, the dimension of the double-trace counterterms would then 
be consistent, curiously, with the existence of a superspace formalism for these terms that 
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FIG. 24: Integrals (40)-(47) appearing in the four-loop amplitude. The notation follows that of 
fig- 1 



preserves manifestly all sixteen supercharges. 



Here we provide details of the cancellations observed in refs. 32| and affirm the double- 
trace power-count bound ( 15.21) for three and four loops. Moreover, by direct integration of 
the results we show that at three loops the double-trace bound is saturated. At four loops, 
we have not evaluated the required integrals in D = 6, so it is possible (though perhaps 
unlikely) that the behavior of the double-trace terms is better than this bound. 
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FIG. 25: Integrals (48)-(49) appearing in the four-loop amplitude. Integral (50) is required for 
the ansatz to match the color-stripped cut at the level of the integrand. However, it does not 
contribute to the color-dressed amplitude. The notation follows that of fig. [20j 
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A. One-loop ultraviolet divergence 



Before turning to the four-loop amplitude, it is useful to review and expand on the lower- 
loop results. Let us consider first the color-dressed one- loop amplitude in eq. ( 12. 5p . For 
each of the contributing integrals it is necessary to specify the corresponding color factor. 
The trace basis, in which the color factors are expressed in terms of traces of products of 
group generators T ai in the fundamental representation, is convenient for discussing the UV 
divergences for an SU(N C ) gauge group. The basis elements are: 



Tr ijk i = Tr (T a! T aj T afe T 



(5.3) 



(In general one may also add to the basis Tr(T ai ), Tr (T ai T a :> T ak ) , etc. For four-point am- 
plitudes the traces of length one and three must appear together; for SU(N C ) they are not 
necessary, because Tr(T a ) = 0.) 

The color factor of the box integral, given in eq. (12. 7p . may be expressed in the trace 
basis as 



Ci B 234 = N c ( Tr 1234 + Tr 1432 ) + 2 ( Tr 12 Tr 34 + Tr 14 Tr 23 + Tr 13 Tr 24 ) . (5.4) 

The lowest dimension in which the one-loop box integral in eq. ( 12. 8 p develops an ultraviolet 
divergence is D = 8. Near D = 8, we have 

pole O [4:71)* e 

From eqs. ( 12. 5 p and ( 15. 4 p it then follows that the divergence of the one-loop amplitude in 
the critical dimension is given by 



4 1} (1,2,3,4) 



SU(N C ) 
pole 



6(47r) 4 e 



N c ( Tr L324 + Tri423 + Tri 243 + Tr 



■1342 



Tr 



1234 



Tr 



1432 J 



+ 6 ( Tr i2 Tr 34 + Tr i4 Tr 23 + Tr i3 Tr 



21, 



(5.6) 



At this loop order, the coefficients of the double-trace terms, relat ive t o the single-trace 
ones, are fixed by a U(l) decoupling identity, or dual Ward identity 10 1| . 

For a general gauge group G the UV divergence is expressed in terms of three independent 
color tensors: the two tree- level tensors and the irreducible one- loop tensor Cf 234 in eq. ( 12. 7\i : 



4^(1, 2, 3, 4) 



G 



pole 



4 1^" 1 



(5.7) 
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where Ca is the quadratic Casimir in the adjoint representation, normalized as in eq. ( ]Bip 
of appendix [Bl Note that the Bose symmetry of the divergence is not manifest in this form. 
As discussed in appendix [HI any one-loop four-point quantity can be expressed in terms of 
these three color tensors. Thus at one loop, in the critical dimension D c = 8, the coefficients 
of all listed color structures diverge and there are no additional hidden cancellations. 



B. Two-loop ultraviolet divergence 



A similar analysis may be carried out for the two-loop four-point M = 4 sYM amplitude, 
which is given in eq. (12. lip in terms of planar and non-planar double-box integrals 1^ and 
j(np) These integrals first diverge in D = 7 and their poles in D = 7 — 2e are 161 ] . 

j(P),D=7-2e| 



y(NP) = 



I pole 



7T 



Sl2 

j(NP),D=7-2e 



20 (4vr) 7 e 

pole 



(5.8) 
(5.9) 



si2 30 (4rr) 7 e' 

For an SU (N c ) gauge group, the planar and non-planar two-loop color tensors C{ 234 and 
^1234^ defined in eq. (12. 131) become, in the trace basis, 



C 



(P) 

1234 



(AT C 2 + 2) ( Tr 1234 + Tr 1432 ) + 2 ( Tr 1243 + Tr 1342 ) - 4 ( Tr 1423 + Tr 1324 ) 



+ 6iV c Tr 12 Tr 34 , 

Ci234' ) = 2 ( Tr 1234 + Tr 1432 ) + 2 ( Tr 1243 + Tr 1342 ) — 4 ( Tr 1423 + Tr 1324 ) 
+ 2N C (2 Tr 12 Tr 34 - Tr 13 Tr 24 - Tr 14 Tr 23 ) . 



(5.10) 



(5.11) 



The full amplitude was originally presented in the trace basis 15|. In terms of and 
1/( NP ), the UV divergence of the amplitude (12. lip is, 

SU(N C ) 



4 2) (1,2,3,4) 



-g b lC 



N 2 C V^ + 12(V^ + V^)) (5.12) 
x ( s i2 ( Tr 1324 + Tr 1423 ) + s 23 ( Tr 1243 + Tr 1342 ) + s 13 ( Tr 1234 + Tr 1432 
- 12 N c (V^ + V^) (s 12 Tr 12 Tr 34 + s 23 Tr 14 Tr 23 + s 13 Tr 13 Tr : 



Inserting the planar and non-planar integral poles (15.81) and (15. 9p into eq. (I5.12p . the UV 
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divergence becomes, 

ro\ su(n c ) n e n JC r / 

4 2) (1,2,3,4) = I (iV c 2 + 20) U 2 (Tr 1324 + Tr 1423 ) (5.13) 

Dole ZU 47T ' e L V 



+ S23 ( Tr 12 43 + Tr 1342 ) + Si 3 ( Tr 12 34 + Tr 1432 ) j 
20iV c [s 12 Tr i2 Tr 34 + s 23 Tr 14 Tr 23 + s 13 Tr 13 Tr 24 ) 



As was the case at one loop, U(l) decoupling identities at two loops |102| relate the 
coefficients of the single- and double-trace structures. These identities provide an a priori 
justification that the double-trace terms diverge whenever the single-trace terms diverge, at 
this loop order. In slightly more detail, the double-trace coefficient is equal to the negative 
of a subleading-color (N®) single-trace coefficient, plus the sum of all three leading-color 
(N%) single-trace coefficients. (See eqs. (4.48)-(4.50) of ref. [1021 ] .) However, in the case of 
eq. ( 15. 13ft . the relevant leading-color sum vanishes by the identity si 2 + s 23 + s i3 = 0. Thus 
group theory enforces the equality of the two "20" s appearing in eq. ( 15. 13ft . 

We may also analyze the color structure for a general gauge group G. In this case, 
the UV divergence of the two-loop amplitude depends on five independent color tensors, 
which we take to be the tree-level and one-loop tensors that already appear in the one-loop 
divergence (15.71) . plus two new independent (irreducible) two-loop tensors. We take the 
latter to be two independent permutations of the tensor Cj^- (All other tensors are related 
to these by repeated application of the Jacobi identity.) In terms of the pole parts ( 15.8ft and 
( 1 5.9ft of the planar and non-planar double-box integrals, the two-loop divergence is: 



4 2) (1,2,3,4)' G 



-C 2 A V {F) {s 12 f a2aib f ba4ai + s 2 3f aiaab f baaa *) (5.14) 

pole I 4 

+ (y(P) + V^)(3C A Cf 2 ° 3 \s 13 + 2C£i(s 12 - s 13 ) + 2C£i(s 23 - s 13 )) 

The five color tensors appearing in eq. (15.14ft form a basis in the space of two-loop four-point 
color tensors (see appendix |B]) . Each of their coefficients diverges in the critical dimension 
D c = 7, and we see no natural combination of coefficients for which the D = 7 divergence 
cancels. 



C. Three- loop ultraviolet divergence 



Let us proceed with a similar discussion of the UV divergences of the three-loop four- 
point amplitude given in eq. (12. 14|) and fig. HJ In this case the critical dimension is D c = 6. 
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FIG. 26: The vacuum-like graphs describing the three-loop UV divergences in D = 6. The large 
(blue) dots indicate that a propagator appears squared, or doubled. 

In the representation of the three- loop amplitude given in fig. HI integrals (a)-(d) and (h) 
are all finite in the critical dimension. Only integrals (e), (f), (g) and (i) contribute to the 
UV divergence, because they have numerator polynomials that are quadratic in the l oop 
momentu m. To extract the divergence we follow the procedure discussed in refs 



11 



103 



29. 



104J, i.e. we expand in small external momenta and keep only the leading terms. In 



this limit, the contributing integrals reduce to, 

j(e)_^ 5l2 y(A) j /(f)_ >% y(B) ) j( g )^ Sl2 y(B) ) 

I (i) ^(s 12 -s 23 ){v^- 1 -V^), (5.15) 

where and V™ are the vacuum-like graphs in fig. [26j These integrals are both UV and 
infrared divergent. As discussed in ref. 29[, to extract the UV divergence we inject off-shell 
momenta at two vertices, thus removing the infrared divergence. The values of their UV 
divergences are: 

1 

pole 



V (B) 



6(4vr) 9 e 
1 



(5.16) 



(5.17) 



pole 6(4vr) 9 e 

The relevant three-loop color factors C^ e \ and are easy to express in the 
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trace basis, 
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3(Tr 1423 + Tr 1324 )) 
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3 (Tr M23 + Tri 324 )) 



(5.19) 



Tr 



1423 



Tr 



1324; 



(5.20) 



(5.21) 



Using these expressions, as well as the relation (I5.15P between the leading poles of the 
integrals I^ c \ I^ s \ and I® and the vacuum integrals and V^ B \ we find that the 
leading UV divergence of the three-loop amplitude in dimension D = 6 — 2e is 



4 3) (M,3,4) 



SU(N C ) 



pole 



2 g 8 K (nI V {A) + 12 N c (V {A) + 3 V (B) )) (5.22) 
x (si2 ( Tr 1324 + Tr 1423 ) + s 23 ( Tr 1243 + Tr 1342 ) + s 13 ( Tr 1234 + Tr 1432 



Inserting the UV pole parts ( 15. 16ft and ( 15. 17ft of and V^ B ' then gives, 



Af (1,2,3,4) 



SU(N C ) 



9 8 IC 



pole 



3(4yr) 9 e 



(A^ + 36C(3)iV c 



(5.23) 



x (^$12 ( Tri 32 4 + Tri4 23 ) + s 23 ( Tri 2 4 3 + Tri 3 4 2 ) + Si 3 ( Tri 23 4 + Tri4 32 



A remarkable feature of the three-loop UV divergence (15.221) in the critical dimension is 
the absence of double-trace terms. From eqs. (I5.18p - (I5.21I) it is clear that such terms exist 
separately for each integral and cancel only in the complete amplitude. This cancellation 
was first noted 32| as a consequence of the calculation described in this paper. Two rather 



different discussions of this property have been presented recently. One approach is based 
on the pure-spinor fo rm aH Sm , both in string theory in the low-energy limit fl and more 
recently in field theory [12J. The other app roach is based on algebraic non-renormalization 



theorems 



341 ]. following up on earlier work 



105. 
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FIG. 27: The leading contributions with numerators quadratic in loop momentum in fig. [J] can be 
absorbed into the two contact diagrams displayed here, up to relabeling of external legs. Further 
rearrangements push the leading terms into the diagrams of fig. [28] 




■Sl2(Sl3 ~~ s 23 



\/(A) 




y(B) 



FIG. 28: At three loops the leading divergence can be rearranged so that it comes from parent 
graphs which are one-particle reducible, depicted here in the S12 channel. The color factors of 
these graphs have no double-trace contributions. The complete contribution comes from summing 
over all 24 permutations of external legs and multiplication by a symmetry factor of 1/4 to remove 
double counts in both graphs. 

Motivated by these results, especially the structure of the leading contributions in the 
string theory analysis of ref. 33|, we rearrange the leading terms in the small- momentum 
expansion of the integrals in fig. H] such that double-trace terms are manifestly absent. This 
representation is found by noting that all terms quadratic in loop momentum, can be placed 
into contact-term diagrams, of the form in fig. [27J This fact suggests that the leading terms 
can be absorbed into the graphs with external propagators displayed in fig. [281 if we include 
a factor of s 12 to cancel the external propagator. After dressing the graphs with color, and 
a factor of (s 13 — S23) so that the graphs have the proper antisymmetry, it is not difficult to 
verify that with the numerical coefficients in fig. [281 the expression 

4 3) (M,3,4) 



1 Jr^l r v[A) T 



n 

S4 



1234 1 VW "+* L/ 1234 1 V( B ) 



subleading, 



(5.24) 



has the proper leading behavior. Here / 



y(A), ly(B) 



Ivcb) and CYzm, C1234 correspond to the 
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integrals and color factors indicated by the graphs in fig. |28j The dropped subleading terms 
are better behaved in the ultraviolet than the displayed leading terms and are not relevant 
for our discussion. The integrals in eq. (I5.24p are essentially just the two vacuum graphs 
and V^ B \ promoted to one-particle- reducible four-point integrals, and dressed with 
appropriate numerator factors. The color factors C^ 34 and are given by 



(14(112(114 



_ yai (22(15 ja^aaa-r jaeaioas j-asanaQ jagaucn jazay A aio j"ai30i4dll J 

= N C (N 2 C + 12)(Tr 1234 + Tr 1432 - Tr 1243 - Tr 1342 ) , 

£,y( B ) _ Ja 1 a 2 a 5 ja 5 aca 7 ja 6 ai 3 aio j?a 3 a 10 a s ja 8 ana 9 Jayja^an j-a 7 aga 12 ja^a^a-is. 

= 12iV c ^Tri 23 4 + Tri 432 - Tr 1243 - Tr 1342 ) , (5.25) 

exposing the manifest absence of double-trace terms in the leading divergence. Because the 
color factors are one-particle reducible, they are proportional to the tree-level color factors; 
they are given by 

Cffi = CV(A) f a ^ b f ba ^ , CXSl = C vC B) f^fba 3 a 4 (5>26) 

These equations may also be taken as the definition of the three-loop scalar invariants C V (A) 
and C V (B). 

The graphs in fig. [281 dressed with the two types of one-particle-reducible color structure, 
reproduce the original representation of the divergence. Thus they account for the absence 
of double-trace terms for the terms with two powers of loop momentum in the numerator. 
Double-trace terms may only arise from subleading terms in the small-external-momentum 
expansion of the integrals in fig. HJ In section IVIBI we will show that the three- loop double- 
trace terms develop UV divergences starting at D = 20/3; therefore the double-trace bound 
suggested in eq. (15 .2p is saturated at L — 3. We note that the leading divergences at one and 
two loops are not associated with powers of loop momentum in the numerator of the parent 
integral (which can be considered contact terms in the UV limit). Thus a rearrangement 
of the leading-divergence contributions similar to eq. (I5.26P is not possible at one and two 
loops. Therefore the single- and double-trace terms have the same critical dimension at 
these loop orders. 

We now turn to the independent question of the divergence structure for a general gauge 
groups G, starting from the representation in eq. (I2.14p . Following a similar procedure as 
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at one and two loops we write the three-loop four-point color tensors in terms of six inde- 
pendent tensors, which we may take to be the five lower-loop ones used in eq. (I5.14p plus 
one irreducible three-loop tensor (see appendix [B]) . Summing over the various color permu- 
tations, the divergence of the four-particle three-loop amplitude in D = 6 — 2e dimensions 
is 

A?\l, 2, 3, 4) G =g 8 JC V (3) (flu f a ^ b f ba ^ + s 23 f^bfba 3 a 4 \ _ ( 5>27 ) 
pole V / 

(Alternatively, one can use eq. (I5.24p to arrive at eq. (I5.27p .) The color tensor of the 
three-loop divergence is proportional to the tree-level color tensor. Only two out of the 
six independent color tensors for a general gauge group G are present in the divergence. 
However, the scalar coefficient involves additional group invariants C V (A) and C V (b), 
defined in eq. ( I5.26p . which do not appear below three loops, 

y( 3 ) = -2 (CV(A)V (A) + 3CV( B )V (B) ) . (5.28) 

They are nontrivial group invariants constructed out of structure constants, with the index 
contraction following the topology of the vacuum diagrams \A A ' and V^ B \- their explicit 
definitions are given in eq. fIBip . and their values for SU(N C ) are provided in eq. (|B3]). Al- 
though they are not reducible to Ca alone, they are related to the standard group invariants 
C A and (d A bcd ) 2 by, 

C 3 

Cy(A) — Cy(B) = — , (5.29) 

o 

Y 2 ^abcd^abcd 

gCVtA) + -CVcb) = -^J- , (5.30) 

where d A bcd is the totally symmetric rank four tensor in the adjoint representation 5 , and 
N A is the dimension of the adjoint representation. However, the invariants C V (A) and C V (b) 
are more natural in our context, because they correspond directly to the color factors of 
vacuum-like diagrams. 

We can make a few observations related to the one-particle-reducible form of the UV- 
divergent terms displayed in fig. [28j We note that the group invariants in eq. ( I5.28P can be 
promoted naturally to rank-three invariant tensors, as suggested by the marked points on 



5 The invariant tensor d a ^ a2a2ai is denned as d A ia2Cl2ai = ± E we s 4 TT A (T a ^T a ^T a ^T a ^) , where 5 4 is 
the set of permutations of four objects. 
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the internal lines of the vacuum diagrams and in fig. |2S1 which denote doubled 
propagators. This connotation suggests that 

l^iabc /~i rabc z^iabc /~i rabc /r ni \ 

<-V(A) = Oy(A) J , O y(B) =Gy(B)J , (^O.dlJ 

are the most primitive yet nontrivial three-loop rank-three color tensors that can be con- 
structed only out of structure constants. All other three-loop rank-three tensors built from 
a single string of structure constants reduce to a multiple of C\f abc . It is remarkable that 
the UV divergence depends only on the three-loop invariants in eq. (15.281) . Moreover, they 
follow a simple pattern: the color invariants and vacuum integrals are in one-to-one corre- 
spondence, with no mixing between the terms, and their relative numerical coefficient can 
be interpreted (via eq. (15.241) ) as having a combinatorial origin. We shall see below that the 
four-loop divergences follow a similar pattern. 

Finally, we comment on the degree of transcendentality 6 of the divergences. In four di- 
mensions, infrared-divergent terms of M = 4 sYM amplitudes, expanded in Laurent series 
around D = 4, exhibit a uniform degree of transcendentality through at least three an d fou r 
loops in the planar case {2], 1^], and through two loops for the full color dependence 106 ] . 



This property is related to the unifor m deg ree of transcendentality observed for the anoma- 



lous dimension of twist-two operators |l07j |. Information on UV-divergent terms in D > 4 is 
more limited. Through two loops, the simple structure of the vacuum integrals ( 15. 51) . ( 15.8ft 
and ( 15.91) enforces a uniform degree of transcendentality. As can be seen from eqs. (15. 16ft 
and ( 15.171) . this is no longer the case at three loops: a non- uniform degree of transcenden- 
tality can and does occur (cf. eq. (15.231) ). It is interesting to note, however, that for gauge 
group SU(N C ) the coefficient of each power of N c does have a uniform degree of transcen- 
dentality, due to the appearance of the particular linear combination + 3 in the 
subleading-color terms in eq. (15.22ft . The same combination of vacuum integrals appears in 
the UV pole for the three-loop M = 8 supergravity amplitude in D = 6 — 2e, ensuring that it 



also has a uniform degree of transcendentality 29(]. On the other hand, for a general gauge 
group, eq. f 1 5 . 2 8 1) does not display any interesting behavior with respect to transcendentality. 
The significance and generality of these facts may be clarified further by evaluating the UV 
singular terms of the four-loop amplitude, as we do below. 



6 Ricmann £ values £„ are assigned degree of transcendentality n, logarithms are assigned degree 1, poly- 
logarithms Li„ degree n and rational numbers are assigned degree 0. 
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D. Four-loop ultraviolet divergence 



The four-loop four-point amplitude is given in eq. (14. ip and figs. 12014251 According to 
eq. (15. ip . we expect the four-loop planar amplitudes to start diverging in the critical dimen- 
sion D c = 11/2. Indeed, an inspection of the parent integrals in figs. [2lH42"5l reveals that many 
of them have numerator polynomials iVj that are quartic in the loop momentum, leading to 
logarithmic divergences in D c — 11/2. For G = SU(N C ) we have no reason to expect color 
single-trace terms to exhibit further cancellations compared to eq. (15. ip . However, because 
all divergences in the critical dimension can be interpreted as arising from contact terms, 
we expect that, as for the three-loop case, divergences from double-trace terms may cancel, 
increasing the dimension in which such divergences first appear; below we show in detail 
how this occurs. 

Of the 49 non- vanishing integrals in eq. (14. ip . 29 diverge in D c = 11/2. Following a 
similar analysis as at three loops, their leading divergences may be expressed in terms of 11 
vacuum integrals, V\ through Vn, shown in fig. [29) 



ii4 -»■ s i2 Vi , lis -»■ -S12V1 , he ->> -s 12 V 2 , In -s 12 V 2 , 

hs — * s 12^1 , 1 19 ~~ >* ~ s 12^2 , ho S\ 2 V 2 , I 2 1 S12V3 , 
h2 S12V5 , J23 ~~ >* S12V5 , I 2 5 — S12V4 , ^26 ~~ ► ~ s 12^2 , 

^29 -> Si 2 (V2 + V4) , I 3 1 -> -S23V2 , 

^32 -> -Sl 2 (Vi - V5) - S23V2 , ^33 -> Sl 2 (V 2 - V 5 - Vq) - S 23 Vx , 

hi -> (s 23 - si 3 )(y 4 - Ve) , J 3 5 -»■ si2^ 3 - s 23 (V 2 + V 3 - V 4 - V 5 + V7) , 

^36^Si2V" 7 , 137^-523(^5 + 2^), J 38 ->• 2S13V2 , 

J 39 -)• 2s 12 (y 2 - V 3 + V 5 ) - S23 (y 2 - 2V 3 ) , /41 ->• s 12 V 8 , 

/42-^Si 2 V 8 , ^43 -S12V5 , / 4 5 -> S 12 (3V 8 - 2V W ) , 



and all other integrals are finite in the critical dimension. 

The color factors of the 50 parent integrals in figs. [2014251 are collected in appendix For 
gauge group SU(N C ) it is straightforward to evaluate these color factors in the trace basis. 



he -> 2^23 - s 13 )(V i -V 5 - 2V 6 + V 9 ) , 

hs -> -s 13 (y 4 + ^5) - s 23 (2V 6 -V 8 + 3V 9 - 2V U ) - s 13 (2V 8 

/49 s 13 (v 3 - X -V A - l -V, + V? + V 8 - ^V 9 - V w + V n ) , 



V 9 ) 



(5.32) 
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FIG. 29: Four-loop vacuum diagrams describing the UV divergences of individual graphs in the 
four-loop N = 4 sYM amplitude. A dot indicates that a given propagator appears squared in 
the integral. For Vio and Vn, the factor of s^q indicates the insertion of s$q = (Z5 + Iq) 2 into the 
numerator of the integral, where lines 5 and 6 are marked in the figure. 

We refrain from including these expressions directly, due to their length. (In appendix |B] 
we decompose the in a basis of color tensors for a general gauge group, and list the 
basis elements in trace basis.) However, the structure of the UV divergence is substantially 
simpler. Upon using the explicit color factors and the reduction to vacuum integrals (15.321) . 
we find that the UV divergence of the color-dressed amplitude depends only on the three 
integrals V, V2 and V^: 



4 4) (1,2,3,4) 



SU(N C ) 



pole 



-6 g w K N 2 C [N 2 C Vi + 12 (V + 2 V 2 + V 8 )) (5.33) 
x (si2 ( Tri324 + Tri 423 ) + s 23 ( Tr i2 43 + Tri 342 ) + s 13 ( Tri 23 4 + Tri 432 



Thus, we find that double-trace terms are absent from the divergence in the critical dimension 
D c , as was the case at three loops. Another interesting feature is that terms independent of 
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FIG. 30: At four loops the leading divergence can be rearranged so that it comes from graphs 
with two spurious external propagators. The graphs with I/S12 propagators are displayed here. 
This rearrangement makes manifest that the color factors of these graphs have no double-trace 
contributions. The complete contribution comes from summing over all 24 permutations of external 
legs and multiplication by a symmetry factor of 1/16, 1/8 and 1/16, respectively, to remove double 
counts in both graphs. 

N c are also absent from the divergence. Only the leading and next-to-leading powers of iV c 
are present in the single-trace divergence. Finally, using eqs. (I5.12p . ( 15 .22 j) and f)5.33p . the 
divergences for L = 2,3,4 all have the form, 

„4 4 L) (1,2,3,4) SU(Nc) = (L - 1)! g 2L+2 fCN^( (l + V planar + 41 V non " planar J 

pole V v N c J N c J 



x (^ s i2 ( Tr 1324 + Tr 142 3) + S23 ( Tr 12 43 + Tr 1342 ) + s 13 ( Tr 123 4 + Tr 1432 
+ ^l,2 x (double-trace-terms) , (5.34) 

where \/ planar and \/ non_planar come from planar and non-planar four-point integrals, respec- 
tively. 

For a general gauge group the leading UV divergence at four loops has a similarly simple 
structure, proportional to the tree-level color tensor: 

G 

pole 



Af\l, 2, 3, 4) G = g 10 K V {4) (si 2 f a ^ b f ba ^ + s 23 faia 2 bJba 3 aA ^ (5 35) 

pole V / 



where 



= 3 (CvM + 2C V2 V 2 + C V8 V 8 ) . 



(5.36) 
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The coefficients Cy 1} Cy 2 and Cy 8 are the group invariants associated with the correspond- 
ing vacuum diagrams. Their expressions in terms of structure constants are collected in ap- 
pendix [Bl This structure is similar to that of the two- and three-loop UV poles in eqs. (I5.27P 
and (15.281) . As at three loops, the four-loop group invariants are not independent; rather, 
they satisfy the following relations: 

C Vl — Cy 2 = — , (5.37) 

I 9 J abed J abed 

Cv 8 — C V2 . (5.39) 

As with the three-loop case discussed above, it is possible to rearrange the UV-divergent 
contributions at four loops into one-particle-reducible parent graphs. This form manifestly 
exhibits the absence of double-trace terms. In this case our representation exhibits two 
propagators in the same momentum channel, which are canceled by numerator factors, as 
depicted in fig. [30j The divergent part of the amplitude then has the simple form 

4 4) (1,2,3,4) = V^El^V^ + I^V^ + ^Vvs] + subleading. 

(5.40) 

where the integrals correspond to the three graphs in fig. (30j and their color factors are 
proportional to the tree-level color factors, 

CJk = C Vi r^ 6 / 60304 , i = 1, 2, 8. (5.41) 

The next step is to evaluate the UV poles for the four-loop vacuum integrals Vi, V2 
and V$ in their crit ical dimension D = 11/2 — 2e . To this end we use the same infrared 



rearrangement 103| (related to the R* operation 1081 ]) that was used 29j to evaluate the 
three-loop vacuum integrals and V^: We inject and remove momentum with 
k 2 ^ 0, at two of the vertices of the vacuum integral, thus transforming it into a four-loop 
two-point integral, which possesses the same UV poles, but no infrared divergences. (The 
infrared divergences arise in the small- momentum limit from doubled internal propagators.) 
We always take the two vertices in question to be connected by a single propagator. Then 
the four-loop two-point integral factorizes into the product of a finite three-loop two-point 
integral and a UV-divergent one-loop two-point integral. 
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FIG. 31: A two-loop integral with central propagator raised to the power n. If n is not an integer, 
G^ 2 \n,D) cannot be reduced to one-loop integrals. 

Some of the three-loop two-point integrals can be evaluated through a similar procedure, 
by factorizing them into a product of two-loop and one-loop two-point integrals. A few of 
the resulting two-loop integrals, such as tho se sh own in fig. [3U are not factorizable. To 



evaluate them we employ the gluing relations 



109| , which require consistency of the various 



ways of factorizing a higher-loop UV-divergent integral into products of lower-loop integrals. 

For example, the diagram V\ has four inequivalent propagators (not related by symmetry), 
leading to four inequivalent factorizations; they are shown in fig. |32j As usual, a dot indicates 
that a given propagator appears squared in the integral. Similarly, the numbers (9 — 3D/2) 
and (10 — 3-D/2) indicate the power to which that propagator is raised, which is determined 
by dimensional analysis of the three-loop integral. All four factorizations sho uld give the 



109 [ 



same answer; this consistency condition is an example of a gluing relation 

The one-loop bubble integral is simple to evaluate. For an arbitrary dimension D and 
powers n\ and n 2 of the two propagators, it is given by jllO | 

d D p 1 (_i)«i+«2 



^-bubble 



where 



ni,n 2 J 



(2tt) d {{p + kf) ni {p 



-G(n 1 ,n 2 )(-A; 2 )- (ni+n2 - D/2) 



(5.42) 



G(n l7 n 2 ) 



T(-D/2 + ni + n 2 )T{D/2 - m)T(D/2 - n 2 ) 
V{ ni )T{n 2 )Y{D - m - n 2 ) 



(5.43) 
(| + 3e,2) 

and (n 1 ,n 2 ) = (f + 3e, 1). Inserting these values into eq. (I5.42p we find that both integrals 



In D = 11/2 — 2e dimensions, for the cases required in fig. [32] we have (ni, n 2 ) 



have the same UV pole, 



C(| + 3e,2) 
G(J + 3e,l) 



1 1 



2ir(j)c 



i i 



2ir(4) e 



O(l) 
O(l) 



(5.44) 
(5.45) 



The finite three-loop two-point integral s ca n be reduced to a set of master integrals using 
the method of integration by parts (IBP) 109J, in particular the algorithm MINCER, which 
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FIG. 32: The UV poles in the vacuum diagram V\ can be determined from the product of a finite 
three-loop two-point integral with a UV-divergent one-loop two-point integral. There are four 
inequivalent ways of doing the reduction, corresponding to different propagators connecting the 
two points at which the momentum is injected. 



is available in FORM 



Hl| . For the two planar topologies occurring in fig. [321 the so- 



called Benz and ladder topologies, the IBP reduction procedure results in integrals that 
factorize into the product of two-loop and one-loop two-point integrals. Equation (I5.42p can 
be applied to the latter. The two-loop integrals can also be reduced to products of one-loop 
integrals, except for the integrals G^(n, D) shown in fig. EU in which the power n to which 
the central propagator is raised is not an integer. However, a gluing relation can be used to 
solve for such integrals in D = 11/2. 

For example, IBP reduction of the top three lines in fig. [32] results in the following relations 
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for the 1/e pole terms, 
6272 



Vi 



-T 5 (- 

25 U ' 



256 



+ 



-4' 2 



2 r 5 (f)-^r 4 (f)rmm; 



i r 2 (f)r(i) 



25 
12352 



^5/3- 



5 
288 



r(|) 



X 



25 V4 ' 5 
G(J + 3e,l) 



r 4 (|)r(|)r(i)-5 UJ u 



r 



l 'G< 2 >(|,£) 



G(f + 3e,2) 

(47T) 11 

G(| + 3e,l) 



(4tt) 



ii 



(5.46) 
(5.47) 



6G ( 2)( 5 n 



^4> 2 



(47T) 



11 



(5.4? 



Equating the three forms for V\ at order 1/e yields 



G (2) (f,f) 



64 
25 



-^r 2 (|)r 2 (|) + ' 



928 T 3 (f)r(i) 
125 



and 



Vi 



■f r 2 (|)r 2 (i) + 



896 r 3 (|)r(|; 



p/l/3^ 

5 U ' 



^r 3 (i)r(i)r ( i; 



+ o(i). 



(5.49) 
(5.50) 

(5.51) 



(47r) 11 e 

The fourth line of fig. 1321 provides a redundant equation. 

A similar strategy leads to an analytic evaluation of V2. As depicted in fig. [331 there are 
four inequivalent ways of factorizing this integral into a product of three-loop and one-loop 
two-point integrals. The same Benz and ladder topologies appear as for Vi, but with different 
configurations of double propagators. After using eqs. (I5.49j) and f)5.50p for Cr( 2 )(|, y) and 
C (|; if)' they all give the same result, 



V 2 



(An) 



11 



^r 4 (i) + ^r3(i)r(i)r(i; 



+ o{i) 



(5.52) 



There are four inequivalent ways of factorizing the non-planar vacuum integral Vs into 
a product of three-loop and one-loop two-point integrals, depicted in fig. [3H This time 
the non-planar three-loop two-point topology is obtained. The IBP equations can be used 
to reduce any integral with this topology down to the master integral "NO m " of the same 
topology, in which all propagators appear undoubl ed. In D = 4 — 2e, gluing relations allow 
one to solve for the value of this master integral 109| as e — > (which is proportional to 
However, in D — 11/2, we find that the gluing relations do not give new information. 
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FIG. 33: Four inequivalent ways of reducing vacuum diagram V2 to a product of a finite three-loop 
two-point integral with a UV-divergent one-loop two-point integral. 



That is, all four ways of factorizing V§ lead to the same expression, 



4 1 V, 



fin 



(47T) 11 2ir(|) 6 



+ 0(1) 



where 



fin - r5( ' } + If r4 ® r ® r ® + 2 N0 ™ • 



(5.53) 



(5.54) 



Although it is not needed for the four-loop M = 4 sYM amplitude, a similar factorization 
and reduction procedure for V9 gives 

1 4 1 V 9 Rn 



(4vr) 11 2ir(|) 6 



0(1) 



where 



vt = r 5 (f) + ^ r 4 (l)r(i)r(i) - 2 NO m . 



(5.55) 



(5.56) 



E. Gegenbauer sums for non-planar three-loop integral 

Although we could not obtain an analytical value for Vg, or equivalently for Vg or NO m , 
we could obtain the result to 13 digits using the Gegenbauer polynomial x-space technique 
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FIG. 34: Four inequivalent ways of reducing vacuum diagram Vs to a product of a finite non-planar 
three-loop two-point integral with a UV-divergent one-loop two-point integral. 



(GPXT) 112[| . This method is based on the observation that, in position space, all prop- 
agators depend only on the coordinates of the vertices they connect and thus, after Wick 
rotation, they may be identified with the generation function of the Gegenbauer (or ultra- 
spherical) polynomials. Expanding them and evaluating the integrals using properties of 
these polynomials reduces Feynman i nteg rals to finitely many (nested) sums. A nice exposi- 
tion of this technique is given in ref . 113| ; in particular the integral NO m in any dimension 
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D is reduced to a triple sum. Using this approach for D = -1 we obtain 



m ~ 50 r(|) k\ 



K = 



, r(n + K + Ji) 1=o 



E 



/! 



(n-Z)! 



(5.57) 



(Z + K + f)(n-/ + K + ~) 



(n + i)(n + k + | )(n + k + f)(n + k + 2) 



(n + 3)(n + k + § )(n + « + f )(n + re + 4) 



+ 



ii 



3n + 4k 



17 



(Z + k + l)(n - / + k + 1) \ (n + 3)(n + k + §)(n + « + f )(n + 2k + |) 

3n + 4k + 6 



(n + §)(n + k + f ){n + k + 2)(n + 2k + 



+ 



n + 2k + ^ 



3n + 4k + f 



(Z + k + §)(n - / + k + § ) y (n + 3)(n + k + f ){n + k + 4)(n + 2k + \) 

3n + 4k + 9 
(n + \){n + k + §)(n + k + 2)(n + 2k + f) 

The sum over I can be done in terms of hypergeometric functions: 
NO m = §^£I>(K,n) 

1 ^ n=0 K =0 L 

(d(k, n) -D(K-|,n+|)- D(k - f , ra) + D{k -2,n+ §))#(«, n) 

+ (-D(«,7l) - £>(k - | j7 1+ I) - £>(k + |,7l) + D(K - 1,71+ §))#(« 



2 ' 
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(5.5? 



where 



v4(k, n) 



r(K + i)r(n + |)r(n + K + |) 



_D(k, 7i) 



k + 



3^2 



n+i 



in 



k! n! r(n + k 



i) + (« + §) 



ii 



(n + |) + 2(K + 



and 



H(k, n) 



K + 1 



3 F 2 (|, k + 1, -n; k + 2, -n - |; 1) . 



(5.59) 
(5.60) 

(5.61) 
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Using this representation and truncating the k and n sums at a value N up to 6500, we 
found sequences of truncated values NO m (iV), such as 

NO m (6000) = -6.197074209444889, NO m (6100) = -6.197095923684655 , 
NO m (6200) = -6.197116937698505, NO m (6300) = -6.197137284819593 , (5.62) 
NO m (6400) = -6.197156996298421, NO m (6500) = -6.197176101462998 , 

which we then fit to a polynomial in 1/N, obtaining 

NO m = -6.198399226750(2), (5.63) 

where the number in parentheses indicates the uncertainty in the last digit. 7 

We also applied GPXT to V$ n and Vg, obtaining similar sums, but with somewhat 
more complicated summands. The sequence for V^ fin converges the fastest with N, but its 
numerical evaluation takes longer. We obtain, 

V 8 fin = 1.428452926283(3), (5.64) 
Kf n = 2.472370645275(3). (5.65) 

To the given accuracy, these values are compatible with eq. (I5.63P and the two analytic 
relations, eqs. (I5.54p and (I5.56p . 

The combination appearing in the subleading- color part of eq. f!5.33j) is 

y r OT , 1 4 1 6.161859216543(3) , c 

* + ^ + v* = (4^21 rgy 1 U • ( 5 - 66 ) 

Thus the ratio of the subleading-color term to the leading-color term in eq. (I5.33j) is 

12 (Vy + 2V 2 + Vg) 44.40538395605(2) 

NfVi = N! • (5 - 6?) 

Amusingly, the large- N c approximation is strikingly bad, if we take the gauge group to be 
SU(3) as in QCD. It is also strikingly bad at three loops, where the ratio analogous to 
eq. (I5.67p can be extracted from eq. ( 15. 23ft and is very similar in magnitude, 36 £3/ 'N% = 
(43.274... )/N 2 c . 

The fact that these numerical ratios are irrational (in the four-loop case, apparently 
irrational) precludes the overall UV divergence from canceling for any gauge group G at 
three or four loops, because ratios of group invariants are always rational. 



7 Recently a m uch more accurate numerical value for t his in tegral has been obtained by Lee, Smirnov and 
Smirnov 



11 4j , using methods similar to those in ref. [115 1 
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VI. UV DIVERGENCES OF SUBLEADING-COLOR STRUCTURES 



As discussed in the previous section, the color double-trace terms are better behaved in 
the ultraviolet than the single-trace terms. Instead of the finiteness bound ( 15. ip . for three 
and four loops the double-trace terms satisfy the bound (15.21) . The double-trace terms of the 
three- and four-loop amplitudes are finite in the dimensions where the corresponding single- 
trace amplitudes first diverge, respectively in D = 6 and D = 11/2. This result implies that 
two-derivative double-trace operators of the form Tr(V 2 F 2 ) Tr(F 2 ) are not renormalized in 
these dimensions. 

Suppose an L-loop cubic parent integral has / powers of the loop momentum in the 
numerator factor iVj. Because it has 3L + 1 propagators in the denominator, it behaves in 
the UV as, 



By dimensional analysis, the number of powers m of external momenta sitting in front of 
this integral (including the four powers in the prefactor K defined in eq. ( 12.61) ) is related to 
/ by m = 2L + 2 — I. The critical dimension in which UV divergences can first appear in 
eq. (16. ip is found by setting DL + / — 6L — 2 = 0. Eliminating / in factor of m, the critical 
dimension is 



By Lorentz invariance, only even powers of momenta give nonvanishing results, implying 
that m is effectively always even. Furthermore, the three- and four-loop integrals have 
at least six powers of external momenta that are manifest in the prefactors, counting the 
four powers in JC and two powers in the "worst behaved" integral numerators Nf, therefore 
m > 6. We showed in the previous section that the divergences in the single-trace terms 
indeed start at m = 6, with no further hidden cancellations. We also found cancellations in 
the double-traces terms in D = 4 + 6/L. If there are no further cancellations, we expect the 
double-trace divergences to start with m = 8, corresponding to D = 20/3 for the three- loop 
double-trace terms, and implying D = 6 for the four-loop double-trace terms. 

An important question is whether the bound (I5.2p is saturated for the double-trace terms, 
or whether further hidden cancellations remain. To definitively answer this question we must 
directly integrate our expressions, in order to extract the coefficients of the potential double- 
trace UV divergences in D = 20/3 and D = 6 at three and four loops, respectively. Below 




(6.1) 



D(m,L) = 4 + -. 

Lj 



(6.2) 
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we answer this question at three loops. 



A. Extracting UV divergences 

rom multi-loop integrals 



A systematic procedure for obtaining the potential divergences 
is based on differentiating with respect to external momenta 103 



104J. Here we follow a 



related procedure, and expand the amplitudes for small external momenta fcj. Formally 
this is achieved by introducing a single infinitesimal parameter e, giving formal amplitudes 
A(ki) —> A(eki) = e 6 a 6 + e 7 a 7 + e s a$ + . . ., where the expansion starts at e 6 due to the 
manifest m = 6 behavior of each integral. At the integrand level this expansion corresponds 
to, 

I {h, £j) -> I{eki, lj) = e 2 v 6 + e 3 v 7 + e 4 v 8 + ... . (6.3) 

In eq. (I6.3P we have dropped the overall tree factor /C, and the v m correspond to sums of 
vacuum-like integrands, depending only on the loop momenta £j, 

V m = ^2 P™p( k i) Vpi^j) > ( 6 - 4 ) 
P 

where the dependence on the external momenta factorizes into polynomials p mp of degree 
(m — 4). The terms with odd powers of e in eq. (I6.3P may be dropped because they are 
zero by Lorentz invariance. To make sense of the expansion under the integral sign, we 
formally take e\ki\ <C \£j\ (where | • | is the Euclidean norm). This is not true in general, 
because the £j are integrated over all values. Fortunately, this error affects only the UV-finite 
part of the amplitude and is therefore not relevant to our discussion. The UV-divergent 
contributions indeed come from the region \ki\ \£j\ (if the amplitude contains no UV- 
divergent subdiagrams). We thus interpret the integrand expansion (16.31) as a formal series 
where the integrals over v m encode the UV divergence of the integral in dimension D(m, L). 
If sub divergences appear, which happens for the potential color double-trace divergen ce at 



four loops, they must be accounted for; a procedure for doing so may be found in refs. 103 



1041 ] . Here we only evaluate the three-loop case, which has no subdivergences. 

When evaluating terms in the expansion (I6.3P one encounters tensor vacuum integrands 
that are not quite in the factorized form of eq. (16.41) . but contain numerator factors in which 
external and loop momenta are contracted, of the form fcj ■ £j. These integrands are easily 
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converted to the factorized form (I6.4p using the identity, 

-+fT^. (6-5) 

This identity is valid for tensor vacuum integrals with two free space-time indices, which 
must be proportional to the metric tensor by Lorentz invariance. (Similar identities for 
tensors with more free indices are easily constructed, but we will not need them here.) 

The various vacuum integrands Vj are distinguished by their propagator structure and 
possibly numerator factors l± ■ £j, and can be represented as Feynman-like diagrams of 
various topologies. There a re m any hidden relations between the integrals, some of them 



generated by IBP identities 1091 ] . These identities complicate the analysis slightly, because 
the representation of v m in terms of the Vj, as given in eq. (I6.4|) . is not unique. To expose 
the identities we can make use of the invariance of the integrals under reparametrizations, 

/L n L 

n 1 ( ek *> * 3 ) = / n 1 £ j + yi c ip £k p) ' ( 6 - e ) 
j=l J j=l p 

where Cj P are arbitrary numbers. After such a reparametrization the e-expansion of the 
integrands look different, but the total UV pole after integration must be the same. We 
equate the different representations of v m integrals, for a sufficient number of reparametrized 
expansions of the form (I6.6|) . 

/ Vm = J^PmpVp = J2p'mp V P= ■ ■ ■ , ( 6 - 7 ) 

p p 

where V p stands for the integral of V p . The resulting system of linear equations for the poles 
of the Vj can be solved easily. In contrast to IBP identities, the relations we find are between 
the pole parts of vacuum integrals that are relevant to our calculation, and no others. 



B. The logarithmic three-loop divergence in D = 20/3 — 2e 

Consider now the amplitude in D = 20/3, which is the lowest dimension for which a 
potential divergence can appear in the three-loop color double-trace terms, and corresponds 
to the v$ terms in eq. (16. 3p . In the expansion of the nine three-loop integrals we obtain 11 
different scalar vacuum integrals. We will not present them explicitly here. However, we 
note that nine of them have a propagator topology very similar to the nine parent graphs of 
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FIG. 35: The three-loop vacuum diagram appearing in the UV divergence of the three-loop four- 
point M = 4 sYM amplitude in D = 20/3 — 2e. Dots indicate the appearance of squared propagator 
factors in the integral. 

the three- loop amplitude, but where the external momenta is set to zero, fc, — > 0, effectively 
replacing a leg insertion by a two-point vertex, and thereby doubling a propagator. (An 
example is provided by V^ h > in fig. |35j) The two remaining scalar vacuum diagrams both 
have a contact interaction. Remarkably, the consistency of the expansion ( 16. 7ft demands 
that the eleven vacuum integrals are all proportional to each other in D = 20/3. The 
proportionality constants are simple rational numbers. We can therefore express the results 
in terms of the single vacuum-like diagram displayed in fig. [35j We obtain the following 
contributions to the three- loop divergence in D = 20/3 — 2e, 

4 3) (1, 2, 3, 4) (D=20/3) = -g 8 K Tr 1234 Qn 3 c (7s 2 + 7t 2 + 6u 2 ) 

Tri234— pole V5 

5 



-N c (23s 2 + 23t 2 - 126m 2 )) , (6.8) 



^i 3) (l, 2,3,4) (D 2 ° /3) = -//CV w Tr 12 Tr3 4 f7A^ c 2 (58s 2 + 3t 2 + 3M 2 ) 

Tri2 Tr 34 -polc \ 

+ 400(s 2 + t 2 + w 2 )) , (6.9) 

where we have isolated the pieces of the divergence proportional to two characteristic four- 
point color structures of SU(N C ): the single trace Tr 1234 and double trace Tr 12 Tr 34 , respec- 
tively. The divergences of all other single- and double-trace structures are simply related 
to eq. ( 16. 8 ft and eq. (16. 9p by crossing symmetry. This proves that no further cancellations 
remain in the double-trace terms and at least for L = 3, the finiteness formula (I5.2p is sat- 
urated. The scalar vacuum integral appearing in eq. (I6.8P and eq. ( 16. 9 p is shown in fig. [35j 
and the UV divergence is given by, 

v "" , = -(4^^io7 + «' (6 ' 10 > 
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We can similarly determine the D = 20/3 logarithmic divergence for a general gauge 
group G. For clarity, we split up the divergence into two contributions: one containing the 
tree color tensors and one containing the irreducible loop color tensors (the latter was shown 
to be finite in D = 6). The result is, 

4 3) (1,2,3,4) (D=2 ° /3) = -^g s K:V^\bfH200C v( n ) (7s 2 -m 2 + 7u 2 ) 

trec-polc oZ(j L \ 

- 16C 3 A (7s 2 + 6t 2 + 7u 2 )) 
- 4 0) (200CV(B)(14s 2 - 87t 2 - 87u 2 ) 

+ \QC\{Qs 2 + 7t 2 + 7u 2 )\\ , (6.11) 



4 3) (1,2,3,4) ' = -^g 8 ICV {h) 21feS 1) Ci(3s 2 + 3t 2 + 58u 2 ) + 770&f C A {t 2 - u 



(73=20/3) 1 



(l)/n(2/Q„2 , Q+ 2 , r e „,2\ , 7 7n^( 2 )^.^2 „,2> 



loop-pole 24^ 

70&i 2) CU(9s 2 + 20t 2 + 31m 2 ) + 1200&! 3) (s 2 + t 2 + m 2 ) 



(6.12) 



where the basis color tensors b[ L ^ can be found in appendix |B] For the three-loop four-point 
amplitude, any color tensor can be written in terms of the six basis tensors that feature in 
the above formulae. Note that none of the coefficients of the basis elements vanish. The full 
divergence is simply the sum of the tree (16. lip and loop (16. 12ft color-tensor contributions. 
Note also that this division into tree and loop contributions is somewhat arbitrary, and 
certainly not unique (unless the loop contribution vanishes, as in the case of the three-loop 
divergence in D = 6); it depends on the choice of basis for the color tensors. (Specifically, 
choosing different basis tensors b[ at L loops affects the coefficients of the b[ L ^ tensors, as 
well as of lower-loop basis tensors, but not the coefficients of higher- loop basis tensors.) The 
split in eqs. (16. lip and (I6.12p does not even respect the Bose symmetry of the divergence. 

By inspecting the above expressions, it is clear that at three loops the coefficients of 
all independent color tensors, for SU(N C ) as well as for a general group G, develop UV 
divergences in D = 20/3. In other words, the double-trace terms for SU(N C ), and the loop 
color-tensor terms for general G, which are finite in D = 6, exhibit no further cancellations; 
their divergences start at m = 8, or equivalently at D = 20/3. 

Here we will not evaluate the four-loop divergence in D = 6. Its evaluation is more 
involved than that of the three-loop case, not only because of the greater complexity of the 
four-loop integrands, but also because of a three-loop sub divergence. A calculation of the 
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divergences in D = 6 would, of course, answer the question of whether there may be further 
hidden UV cancellations in the four-loop double-trace terms. 



VII. CONCLUSIONS 



In this paper we obtained the four-loop four-point amplitude of M = 4 super- Yang-Mills 
theory, i.e. both the planar and the non-planar contributions, in terms of a set of 50 loop 
integrals (one of which has vanishing color coefficient and also integrates to zero) . The planar 



131 ] . Here we presented the more complicated non 



part of the amplitude was already known 
planar terms. At the core of our computation lies the unitarity method 
the method of maximal cuts 



2n 



13 



and in particular 



64J. For the planar amplitude, hidden symmetries [4J lead 



to major simplifications. No comparable considerations are presently available for the non- 
planar terms. 

An important problem is therefore to develop new tools for constructing multi-loop in- 
tegrands, which are valid also for non-planar amplitudes. In this paper we extended and 
developed graphical rules applicable to such amplitudes, to speed up their construction. 



The simplest of these rules is a generalization of the rung rule [15|, [16], which allows us 
to write down all contributions having two-particle cuts directly from lower-loop results. 
However, not all terms have two-particle cuts. A more powerful tool, which partly bypasses 
this limitation, is the box cut, which yields terms having four-point subgraphs. This ru 
has a simple graphical formulation. Another rule, based on the color-kinematic duality 



allows us to generate many non-planar contributions from much simpler planar ones. In 
fact, it appears that these relations can be applied directly at loop level with no cut condi- 
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Although the 



tions imposed, though this has been confirmed only through three loops 
graphical rules presented in this paper do not determine all contributions, they determine 
most terms, allowing us to focus on the remaining ones, for which we employed generalized 
unitarity, in particular maximal and near-maximal cuts. Nevertheless, identifying further 
tools and structures would be extremely helpful for future studies. 



The results in this paper, as well as those of refs. [10J, [11|, |29|, [64 1, provide a wealth 
of information on the explicit form of planar and non-planar amplitudes in M = 4 sYM. 
However, the structures we are seeking, which could allow the development of new tools, 
are obscured by ambiguities in assigning contact terms to the parent graphs. Each form of 
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the amplitude may expose one property or structure while hiding others. For example, the 
representations of the complete three- and four-loop four-point amplitudes presented in this 
paper do not highlight the fact that double-trace terms are better behaved in the UV than 
single-trace terms. As we discussed, a nontrivial rearrangement is necessary to manifestly 
expose this feature. As another example, nontrivial rearrangements would be needed to 
expose a recently proposed loop-level color-kinematics duality in the four-loop four-point 
M = 4 sYM amplitude. It is quite likely that other important features will be revealed 
through appropriate reorganizations of the amplitudes. 

Generally it is simplest to construct an ansatz for the amplitude in four dimensions, where 
powerful helicity and on-shell superspace methods can be used. Because we are interested 
in the higher- dimensional UV structure of the theory, it is important that our construction 
of the four-loop amplitude be valid in higher dimensions. In this paper, we evaluated a 
complete set of cuts in four dimensions, but only an incomplete set in D dimensions. Strong 
checks that our expressions are valid in D dimensions come from two-particle cuts, box cuts, 
and the color-kinematic duality 36] . The planar contributions to the four-loop four-point 



amplitude have been evaluated in D dimensions, subject only to the mild condition that no 



terms violate the expected power count 



13] . It is nevertheless still important to compare 



our results against a complete evaluation of the D- dimensional cuts. At least in D = 6, 



efficient tools for doing so now exist 



sa-lssi- 



A key reason for computing the four-loop M = 4 sYM amplitude is to study its UV 
behavior as a function of dimension. Our expressions are manifestly finite for D < 11/2 in 



accordance with the expectation 16|, |30| that L-loop M = 4 sYM amplitudes are finite for 



D < 4+6/ L. Direct evaluation of the integrals reveals that the theory is indeed UV divergent 
in D = 11/2 for a general non-abelian gauge group G; we have therefore demonstrated that 
no hidden UV cancellations remain for generic G, through at least four loops. 

In the course of our analysis of the color structure of UV divergences, we have found, 
however, that certain terms are more convergent than the general expectation: in particular, 
for G = SU(N C ) the terms with double-trace color factors have no divergences in the critical 
dimension for which single-trace divergences first appear, at three and four loops {32]. This 
curious improved behavior has been discussed from the vantage point of the pure-spinor 



formalism in string theory 



331 ] and field theory 12 L as well as from the standpoint of field- 



theoretic algebraic non-renormalization theorems 34J. In the present paper, we give a field- 



75 



theoretic clarification, motivated by the structure of the leading contributions in the string 



theory analysis of ref. 



33] : It is possible to rearrange all leading-divergence contributions 



of the amplitude into a one-particle-reducible form, which has propagators depending only 
on external momenta. The color dressing of such terms forbids the presence of double- 
trace color factors. In contrast to the string-based approach, this also clarifies the absence 
of subleading-color double-trace terms. This rearrangement can first be performed at three 
loops because this is the first loop order where inverse propagators can appear in numerators 
of integrals. The first potential divergence of double-trace terms is then at D = 4 + 8/L, 
in contrast to the D = 4 + 6/L divergences for single-trace terms. We showed that the 
D = 4 + 8/L bound is saturated at three loops; it remains an interesting open question 
whether it is saturated at four loops. 

An equally thorough understanding of the UV properties of M = 8 supergravity is absent. 
However, the complete four-loop four-point amplitude constructed in this paper helps shed 
light on this issue. This is another important motivation for our work. The M = 4 sYM 
amplitude is the basic input into the construction of the corresponding M = 8 supergravity 



amplitude 
unitarity 
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bllowing the strategy of refs. 10 



16 



: as a consequence of generalized 



62] and of the KLT 35] and graphical numerator double-copy relations 



cuts of N = 8 supergravity amplitudes can be expressed directly in terms of the cuts of 



M = 4 sYM ami 



UV properties |9Nl 1 



itudes. Although M = 8 supergravity is known to have surprisingly good 



16 



29. 



30 



43 



44 



105] and may even be U V finite l9| , its U V behavior 



beyond four loops is still unclear. Recently, a consensus has formed 33 



34 



46 



47 



116 



11S | 



that supersymmetry alone cannot protect the theory beyond seven loops in D = 4 and 
that an additional mechanism must be at play if the theory is finite. Clearly, further high 
loop studies would be helpful for resolving this issue. Apart from its direct relevance for 
answering the question of ultraviolet finiteness of M = 8 supergravity, M = 4 sYM offers a 
much simpler arena for sharpening our understanding of such divergences in supersymmetric 
theories, as well as for explicitly testing general arguments. 

The results described in this paper may also shed light on the infrared singularities of 
gauge-theory amplitudes, by allowing tests of proposed structures for the soft anomalous- 
dimension matrix j^, [|3l| l^] and the cusp anomalous dimension. The evaluation of the 
non-planar integrals appearing in the expression of the amplitude in D = 4 — 2e is a necessary 
step for the extraction of its infrared divergences. Unfortunately, such four-point integrals 
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are notoriously difficult to evaluate; they have not been computed even at three loops. Thus, 
the use of our results for this purpose must await the development of new techniques for 
evaluating higher-loop non-planar integrals. 

In summary, the results and tools presented here should provide new handles on a num- 
ber of important unanswered questions, including the structure of multi-loop infrared di- 
vergences in gauge theories and the UV properties of gravity theories. They may also help 
expose hitherto unexpected structures in these theories. 
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Appendix A: Sample evaluation of a non-planar cut 

In this appendix we describe the details of the evaluation of a nontrivial cut which appears 
in the construction of eq. (14.11) . The cut presented here can be used to confirm large parts of 
our construction, and is sensitive to the most complicated non-planar terms in the amplitude. 
Seven-particle cuts, which break the four-loop amplitude into four- and five-point MHV and 
MHV tree amplitudes are especially helpful in our construction, because they are sensitive 
to all terms actually present in the amplitudes, and yet are relatively simple to work with. 
As mentioned in the text, we checked that our expression for the amplitudes matches all 
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such cuts. 





(a) 



(b) 



FIG. 36: A nontrivial non-planar cut at four loops. The cuts (a) and (b) represent the two distinct 
internal helicity configurations contributing to this seven-particle cut. The blobs representing MHV 



trees are labeled by a "+" and the blobs representing MHV blobs are labeled by a "— " . For a four- 
point amplitude either assignment can be made; the assignment in the figure makes the supersum 
calculation most transparent. 



Here we present the steps necessary to test eq. ( 14. ip and the numerator factors assigned 

This cut has already been 



to the integrals, using the seven-particle cut displayed in fig. 

discussed in considerable detail in sect. 6C of ref. {z^], so here we will simply quote the 



result, and focus instead on the cut of our result in eq. (14. ip . As noted in fig. [36J the 
contributions to the cut comes in two distinct sectors, depending on the configuration of 



MHV and MHV assignment of the tree amplitudes composing the cut. The fact that all 



factors of tree amplitudes are either MHV or MHV leads to compact expressions for these 
two components. From ref. [3], we have the simplified result: 



JMa) = J2Af HV (l-,k,k,hA + )Af^(-h,k,-l 2 ,h,l 5) 



states 



x Ai{-1 4 , -h, 2", -l 7 ) M~h, ~h, h, 3 + ) 
((12)[/ 7 2][/ 6 3][14]) / 
1 



[41][l/ 6 ][/ 6 /2][/ 2 /i][/i4] 
1 



3Sb) 



(/i k) (kh) (hh) (hk) (kh) [3z 5 ][y 6 ][W7]M [2i 7 \[hk}[hh][h 1 A 



[ MHV { 



-li, k, —I2, h, h) 



states 



X A^(— Z4, — 13, 2 , —If) Ai{— / 5 , — /g, ^7, 3 + ) 
,4 1 



((12) [2 3] [Z3/4] [kh]) 



(41) (lZe) {kh){hh) (Zi 4) 
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1 1 1 



(Al) 



x 



This simple result is determined by the purely gluonic configurations crossing the cuts; in 
the first case there are seven such configurations and in the second eight. 
The above two expressions are related by complex conjugation: 



This relation guarantees from the outset that after dividing by the tree amplitude, the sum 
of the two expressions is a function of only Lorentz dot products, with all spinors or Levi- 
Civita tensors dropping out. (This property is special to four-point amplitudes and does 
not hold for higher-point amplitudes.) 

The cut in eq. f lAip needs to be compared to the cut of our expression for the four-loop 
amplitude in eq. (14. ip . To identify the possible contributions to this cut we write out all 
possible tree graphs with only three-point vertices, for the tree amplitudes composing the 
cut. These are depicted in fig. [36j (We need only track graphs with three- vertices, following 
our organization of contributions according to such graphs.) The graphs labeled by (a) and 
(b) correspond to the two five-point amplitudes in the cut, and those labeled by (c) and (d) 
to the four-point amplitudes. Two graphs were dropped from (a), in which external legs 1 
and 4 fuse together, because they would generate three-point subgraphs, which we know do 
not occur. 

After sewing together the four trees, we can identify the integrals and labelings of internal 
and external lines that contribute to this cut. Of the 60 possible such products, 22 may be 
ignored as they do not appear in the list of integrals composing the answer. The unnecessary 
ones would have triangle or bubble subgraphs, and would violate the known no-triangle 
property of M = 4 sYM at one, two or three loops. Organized in this way, the contributions 
to the cut in fig. [36] are listed in table [J which is structured as follows. The first column 
represents the product of trees - - (ijkl) standing for the product of the zth tree from 
fig- l37T a). the j'th tree from figure 137T b). the fcth tree from fig. l37T c) and the /th tree 
from fig. l37T d). The second column is the number n of the integral X n for which {ijkl) 
represents a cut contribution. A given integral may appear multiple times in the same cut, 
Several such products may represent cuts of the same integral. The third column contains 
the mapping of (internal and external) momenta between figs. I20H251 and the momenta of 




(A2) 
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k 1-- 

h 

h 4+- 



.h 4" 



■/o 



/ 2 



(b) Z 4 
ZiJ- 



ZA 



. Z5 h- 



■ h 



■ h I 




h- 



1 1 1 



(d) L 



h 



h ■ 

■ 3+ ^5 



FIG. 37: Tree structures building up the cut in fig. 

the tree amplitudes assigned in figs. [36] and |371 The last column represents the numerator 
factors from figs. [20H25| evaluated for the cut conditions specified in fig. |36j We specify the 
mapping only for the external momenta or when the numerator contains a given internal 
momentum. For convenience we have introduced the notation, 



kij &i ~t~ kj , li j \{ ~\~ lj , 



h,j,m 



j Zj lj , 
- /j lj l m . 



(A3) 



In general, subscripts separated by commas denote sums of momenta, while a bar over the 
subscripts denotes that the sign of the corresponding momentum is minus. 

In table [J there are a total of 38 contributions to the cut. The first such contribution is 

r (i32i) = K + h~ h) 2 [siz{h + k~ h) 2 ~ s 2 3(h + h + l 7 f] 

(fa + h) 2 (h + fa) 2 (h - /i) 2 (/ 3 + h) 2 (h + h) 2 (h ~ fa) 2 ' 1 ' 

where the numerator is given in table [J and the denominators can be read off from the 

propagators of the diagrams in fig. [371 We have put back the overall prefactor K, defined 

in eq. (12.61) . The reader may check that the sum over all 38 contributions adds up to 

give the sum of the two contributions in eq. ( 1A1I) . It is convenient to make this comparison 

numerically, but it is also straightforward to carry out analytically after imposing momentum 

conservation. 
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TABLE I: The numerator factors and tree structures appearing in the seven particle cut of fig. | 
The first column gives the tree structure according to fig. [36] and the notation described in the 
text. The second column gives the integral number n whose cut gives the tree structure. The 
third column gives the relabelings needed to go from the integral labels to the cut labels. The last 
column gives the numerator factors in the labeling of the cut. The associated propagator factors 
can be read off from the tree structures. 



Trees l n 



Momentum relabeling 



Numerator factor 



(1321) 34 



fcl ->■ fa, k 2 -> fci, 
fa — > ki, ki — >■ k 2 , 
h — > fa — h, h — > h 



h + l 5l ) 2 [si3{k 2 + l 5 j) 2 - s 23 (k li3 + l 7 ) 2 ] 



(1322) 18 



fci ->■ k 2 , k 2 

fa — > fa, ki 

h ->■ k 3 + l 7 , 

h -> fci + h, 
l 7 — > I4 



fa, 
fc 4 , 



S2 3 (fcl,3 + ll?{kl A + kf 



(3412) 18 



fci 
fa 



fa, k 2 
fa, fa 



fci, 
fa, 



su{fa-h) 2 {k 2 -l 7 f 



h- 


> _ ^6) 


Z 6 -> 


-k, 


h- 


>-/ 2 






fcl - 


-> fc2, 


fc 2 


-»• fc3 


fa- 


-)■ fc 4 , 


ki 


— >• fci 


h- 


> is,5, 


k~ 





(1522) 20 



«23 ^,4,5(^1,3 + ^7) 



^7 -> fc3 + ^7 



(3212) 20 



fcl- 


■> fa, 


fc 2 - 


■> fcl, 


fc 3 - 


■> fc2, 


fc 4 - 


■> fc3, 




► - Z l,5> 






z 7 - 


► -ie 







si4 Xb {fa-hf 



81 



Trees X n 


Momentum relabeling 


Numerator factor 


(2212) 25 


k\ — > fc 2 — > ki, 
ki — y ko. kA — y ki. 
^6 - > ~k, k 
k — > —^i,5, 
k ~h 


sulfa ~ kfl\ X5 ~ si 4 (/i, 5 + ^) 2 ] 


(2412) 25 


ki — y ki, ki — > ki, 
k 3 -»■ fc 2 , fc 4 -»■ fc 3 , 
^6 — ^ ~h,E: k — > 

Ig — > —I2, 69 — > —I7 


-s M [(A; 2 - Z 2>6 ) 2 (/ 5 - k a f - sufa - l 2 f] 
+ sull t ifa - kfif 


(1511) 33 


k\ — > ki, k 2 — > ki, 
k 3 ->■ fc 3 , fc 4 -> fc 2 , 

^5 — ^ ^3,5; k — * ^4; 

^9 ~~ ^ ^3 ~~ 
/10 k,7i 
In — )■ ki + k, 

ki ~k 


si^lsfaz - k? + 11 jfa + h^) 2 } 

- I 2 [shqCA-i , - h) 2 - ^i(ki 4- /„ E ^ 2 4- St,/ 2 1 

t 3,4,5[ 6 23^1,3 f 5j a 13 ^3 ^ t 4,5,7j 1^ i 13 t 4,7j 


(1312) 33 


ki — > ki, k 2 — > ki, 
k 3 ->■ fc 3 , fc 4 -> fc 2 , 

Z5 — >■ / 4 , ^6 — ► k,5: 

I7 — y — 17. — y 1^. 

^9 ~~ ^ ^3 + ^7) ^10 ~~ * ^3) 

111 — 7" fci + ^6) 

112 — > ~k 


^3 4 5 [ S 14 (&1,3 + h) 2 ~ Sisfa + Z37) 2 ] 

- si 4 (A; 4 + /4) 2 ^7 
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Trees X n 


Momentum relabeling 


Numerator factor 


(1412) 33 


ki k 2 , k 2 -)■ k 3 , 

^3 — ► &1> ^4 - > ki, 

h ~hi h —> —^i,5) 

— y ki + /fi. Zs — ^ /c:. 
/g — > —Iq, Zio — >■ — Zl, 

111 - ^ — ^7? 

112 k 3 + l 7 


~(h ~h- h) 2 [sn(h + hf - s 23 (k 3 - k) 2 } 
-s 23 (A; 2 -/ 2 ) 2 (A; 1 + / 5 + / 6 ) 2 


(1311) 36 


ki ->■ fc 2 , & 2 -> fc 3 , 

&3 — ► fci, &4 — > ki, 
h — >■ Z4, ^6 ~~ ^ ^5) 

h ^ k 3 — l^, Ig — y lij, 
ho ~~ ^ —hi hi — > h,5i 
hi ~h 


8u(h + h) 2 (k 2 - h) 2 - l 2 5 jSu(ki + li) 2 

+ ^3,5 S 14[ S 13 ~~ (&1 + hfi) 2 ] 
+ (ki 4 + /fi) 2 fSl4(A;i + U a) 2 

+ s 13 ((k 3 + h A7 ) 2 - sii)] 


(1212) 36 


k\ — > ki, k 2 —> ki, 
k 3 ->■ k 3 , h ->■ fc 2 , 

// 7 — r — A 7 . Ac — 7 — L ft . 

/g — > —li, ho ~^ hi 
hi —hi 
l\2 h + l 7 


-(k 2 - h) 2 [su(ki + hj) 2 + si 3 Z? 6 - S14S13] 

+ S U (k 2 + k,s) 2 {h + hfi) 2 

— ^I i 5 s i4[(^i + h,e) 2 ~ | s i3] 


(1512) 36 


ki -)• k 2 , k 2 ->■ k 3 , 
k 3 — > ki, ki — > ki, 

h — > h,5i h ~hi 
h — > ki + /«, 
/ 8 -> k 3 + Z 7 , /g -»■ Z 3 , 
Zio — > hi hi hi 
I12 — > ~h 


^ 3 2 5^i4[si3 + 2(A; 1 + / 6 ) 2 -2/ 5 2 ? ] 
+ s 23 [l 2 5 - 7 (ki + Z 3i5 ) 2 - (fci + Z 6 ) 2 (fc 2 - / 3i5 ) 2 ] 
+ (ki A + h) 2 [s 23 (s 13 - (ki + hj) 2 ) 
- si 3 (k 3 + l 3J ) 2 } 



S3 



Trees X„ 



Momentum relabeling 



Numerator factor 



(1521) 32 



fci -} fa, k 2 ->■ fci, 
fc 3 -> fc 4 , fc 4 -> fc 3 , 
h h, h — > fa — h 



(fa + hj) 2 [si2(k 2 + hj) 2 - s 12 (fa - h 
- s 1A (fa + h) 2 } 



(2411) 21 



ki — > fa, fa — > fa, 
fa ->■ fc 2 , fa ->> fc 3 , 



S14/2 7 (fc 3 — il) 



(3411) 22 



fci fci, fc 2 —> fa, 
fa -> fc 2 , fc 4 ->■ fc 3 , 



Sl4^5 1 2(^3 ~ ^6 



(2211) 23 



fci ^4, fa — > fa, 

fa ->■ fca, ^4 -> fcs, 

Z5 — £7, — > —I2, 

h -> -h 



S14/2 7 (fc 3 — i\Y 



(3211) 23 



fa —7- fci, fc 2 — > k^, 

fa ->■ fc2, &4 -> fc3, 

/ 5 —7- — ^ ~~ ~h, 



S 14^5 1 2 (^3 ~~ ^6 



(1211) 37 



fci ->■ fc 3 , fc 2 -» fc 4 , 

fc 3 ->■ fci, fc 4 fc 2 , 

^5 _ > — ^7? ^6 - > — ^1,5) 
/ 7 -> fc 4 + /l, 

z 8 -> fci + h, 



(&2 + Zl, 5 ) 2 [s34(fcl,4 + Z 6 ) 2 - Si 4 (fc 4 " 

+ (fa ~ lr) 2 [si3(fa,4 + Z1) 2 - Sl 4 (fcl 

+ Si 4 (fc2 + /l, 5 ) 2 (^2-Z7) 2 



- hf] 
+ kf) 2 } 
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Trees X n 


Momentum relabeling 


Numerator factor 


(1411) 37 


ki fci, &2 ->■ fc 3 , 
A; 3 -» fc 4 , fc 4 ->■ fc 2 , 

^5 — ^2 ? ^6 — >■ — ^7; 
^7 - ^ ^5? ^8 ~> — ^1) 
^9 — ► ^1,5 


(fa, _ l 7 ) 2 \ S i*(ki -U) 2 - sul 2 J 

+ (Afe - / 2 ) 2 [si 4 (fc4 + hf ~ S M ll 2 ] 

+ s u (k 2 - h) 2 {h - hf + s 13 sulj 5 


(2522) 13 


fci fci, fc 2 — > k^ 

^5 - ► — ^2,6, 
*6 ^3 + h 


s li(k3 + ^7,2,6) 2 


(2311) 27 


ki — ;> fci, fc 2 — > fc 4 , 

A; 3 -»■ fc 2 , fc 4 -»■ fc 3 , 

^7 ~^ ^5? ^8 — ~h, 
Ig — > —Iq, Zio — >■ — Z 4 , 
^11 — >■ ^7? ^12 - ^ — ^3,5 


Su[sul^ 5 — Sn(k 3 — hfi) 2 } + Su(k 3 — l2fi) 2 [ll : Q — l 3t5 \ 


(2512) 27 


fci — > fci, fc 2 — >■ fc 4 , 

fc 3 ->■ fc 2 , fc 4 ->■ k 3 , 

I k — ¥ — lo a. 1 t — ¥ — / 1 . 

Zg — > —I2, Ig — > —Iq, 
ho ~~ ~h,5i hi ~h: 
l\2 ~h 


4s( k 3 - hfi) 2 - S 23 [s 23 £| 7 + (fc 3 - hfif^fifi] 


(2511) 26 


ki — > ki, k2 —> ki, 
k 3 -»■ k 3 , h -»■ fc 2 , 

fc^ — r — 19 fl« fcfi — r — £fi, 

h ~~ ~hi h ~~ >* h,Ai 
h h,7, ho ~hfi> 
hi h 


sidljKh + h) 2 + (fc 3 + h) 2 - 1\ - 1 2 ] 

— sfil 2 7 — si 4 i| 6 (fc 2 — hfi) 2 + si 4 Z 46 (A; 3 — hfi) 2 
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Trees Z„ 



Momentum relabeling 



Numerator factor 



(2321) 26 



ki — > ki, k 2 — > k 4 , 

k 3 -»■ k 2 , h -»■ fca, 

Z5 —7- — /1, ^6 — >■ — ^4, 

h ~h,5i h ~h, 

h h, ho ~hi 
hi -> £7 



S 14^5 2 ~~ s 14(^2 



3,4,5 



(2521) 26 



k\ — > k±, k 2 —> k±, 
hi ->■ fc 3 , fc 4 ->■ fc 2 , 
£5 — > ~ hi h ~h,5, 
h ~~ ~hi h ~~ ~hi 
h — > k 2 — hi ho — > ~hi 
In — > k 3 — h 



Su(ki — h) ^3,4,5 — s 14 {k 2 — l 2 - 

+ s li (k 2 -i 7 f[(h + i 1 ) 2 -q] 



(2312) 26 



k\ — > k^ 


k 2 ->■ fci, 


h -> h, 


/c 4 -> fe, 


h — * ~hfii 


^6 ~~ >* ~hi 


h ~h,5i 


h ~h,i 


h h, 


ho —hi 


l n -> -l 7 




k\ — > k±, 


k 2 ->■ fc 4 , 



S 14^3,4,l ~~ S 14(^3 - ^2,6) 2 £3,5,6 

+ s M ll 6 (k 2 - hf 



(2322) 12 



k?. ->• ko 



kd — >• ki 



(h + /2,6, 7 ) 2 



h ~h,6i h ~ ► ki + h 



(2122) 40 



fcl 


-> hi 


A; 2 ->■ fc 4 , 


A-3 


-> fcs, 


A; 4 ->■ fc 2 , 


h 


-> £3, 


^6 ~ ^ £1,4 


k x 


-> fcs, 


£; 2 ->■ k 


h 


-»■ fc 4 , 


A; 4 — t- fc 



„2 72 
6 23 t l,3,4 



2- 



(1122) 42 



h^ h + l 7 , 
h ->• fei + ^ 



si 4 (&i, 3 + / 7 ) 2 (^i,4 + /e 
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Trees Z„ 



Momentum relabeling 



Numerator factor 



(3112) 42 



ki — > hi, k 2 —> fc 4 , 
k 3 -»■ k 2 , k 4 -»■ fc 3 , 



Sl4(/C3-^) 2 (^2-/7)' 



(2121) 44 



fcx — > fci, fc 2 — > k 4 , 

A; 3 ->■ fc 2 , fc 4 ->■ fc 3 , 
is -> -ii, 
Z 7 -> fc 2 - i 7 , 



si 4 [(fc 2 - h)\h - hf - (h - h)\k 2 - h, 7 y 



(2112) 44 



ki — > k 4 , k 2 — > fci, 

&3 -> &2, ^4 fc 3 , 
^5 ~~ ^ — ^2,6? h ~^ hi 

h^k 3 + l 7 



—Si4(k 2 — hfif(k 3 — h,6jf + si 4 (A; 3 — l 2) e) 2 l^ 3{ 



(3111) 43 



ki — > k 4 , k 2 — > ki, 
k 3 -»■ fc 3 , fc 4 -»■ k 2 , 

h ~~ ^ ~hi h ~~ ^ ~hi 
h —k\ — h, 
h ~~ h,2 



su(k 3 - hf(k 2 - h 



(2111) 45 



ki — > k±, 
k 3 ->■ fc 2l 
is -> -ii, 

^9 — -^4,7, 
ill — ► i 2 ,3; 
il3 ->• ki + h 



k 2 -)■ fc 4 , 
fc 4 -)> fc 3 , 
i6 - ► — i 2 ,6j 
is - > hi 
ho hi 

il2 — >■ — i7) 



-si 4 p 3 -ii) 2 i| 5 + (A; 2 -/ 1 ) 2 / 



H 2 -1 



2(2 



+s 14 (/c 2 - h) 2 (k 3 - h) 2 + s u (k 2 - hfl 



2.3 



(1121) 46 



k\ — y k 3 , k 2 — > k 2l 
k 3 — > k±, /c 4 — >• ki, 

h — > k\ + hi h ~~ hi 
h — y k 2 — hi h ~hi 
h ~+ ~h 



(h ~ h) 2 [s23(k 4 + h) 2 ~ sia(fci l4 + hf] 
+ (k 2 - h,7) 2 [si 3 (k 1A + hf - s 23 (k h3 + l 7 
- s 13 (h + h)\k 2 - l 7 f 
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Trees X n 


Momentum relabeling 


Numerator factor 




ft/1 7 ft/ 4 , 


ft/ 2 r ft/i, 






r ft/3, 


£• . i. £>„ 

ft/4 7* ft/2, 










Suffci + Z S ) 2 (A;3 - Z1) 2 -s u fe - Z 2 ) 2 (A;4 + L) 2 




/c — \ hi -4- 1-7 

16 t A-3 17, 




+Sl4ffel + Z 3 ) 2 (Zi + Z 6 ) 2 +Si 4 (/C3 - Zi) 2 (/C3 + Z3 7 ) 2 




]i — y — Zo 

fc 7 ' •'z? 


/„ i. 7. 




(1112) 48 




Zm ^ Zo 


-Sl 2 / 2 , 6 (/«3 + ^3,7) 2 +514(^3 + /7) 2 (^3 ~ ^) 2 




hi — > —h,i, 


Z12 — > h,zi 


+si 4 (A;3 + Z 7 ) 2 (/«2 + Z 2 O 2 -513(^3 + Z 7 ) 2 Z 2 fi 




l\3 — > ^i + h, 




-S 14 (/C 1 + Z 6 ) 2 Z? 47 -Si 4 (/€ 3 + ^7) 2 ^,7 




ll4 — > — Z7, 




-s u (k 3 + hfl\ 3 




/15 — )■ fc 4 + /1, 








Zi6 — >■ Z 4j 7 








fci — > fc 4 , 


k 2 -¥ kx, 






7 1 7 

k 3 ->■ fc 3 , 


7 1 7 

k 4 -)• fc 2 , 


Sl4l K 2 f 7j l K 3 f 6j + Sul^^ + 7 3 J \K% — Iq) 




7 .7 I 7 

£5 — 7> /Cl -h fcg, 


7.7 

£ 6 — > — L Q j 


— ZsoSiqfZci -|- Zfj) 2 + ti nS^^ik/i -\- Z1) 2 




7 7 


7 7 

^8 -> -Z7, 


+ Z| 7 si 4 (/ci + Z 6 ) 2 - Z| 7 si 2 (A/4 + /1) 2 


(1111) 49 


f 9 - ' ^3? 

Zn —7- — Zi, 


7 v 7 
tlO —7" 74,7; 

Zl2 ~ >■ ^5) 


- l\ 6^13(^2 - M 2 - Zi,[si3(A;4 - Z 2 ) 2 


+ si 4 (s 13 - (hi + Z 5i6 ) 2 )] - s 14 (/€ 2 - Z 2 ) 2 (A;i + Z 5>6 ) 2 




Z13 —7- — Z4, 


Z14 — > Z 2j 3, 


+ si 3 (/«2-/2) 2 (A;i,4 + ^) 2 + si3(A;3 + ^7) 2 ( fc i,4 + ^) 2 




^15 — >■ 'l,4i 




+ Z 2 7fsi3(fc4 + Z 2 ) 2 -Sl4((Z«1 +Zq fi ) 2 + (Z«4 + Zi a) 2 )! 




/i6 -h- k 


4, 


— {k\ + ^6) 2 [ — 513(^2 — ^1,4) 2 + ^5 S 13 — ^I,4,7 S 12] 




In — ► h — ^3 




+ l\ 5 S13S14 + Z| 3 Z 4 7 s 14 



ss 



Appendix B: Color factors: tensors and Casimirs 



In this paper we encountered several types of scalar color factors, or Casimir invariants, 
as well as group invariants associated with various vacuum graphs. These factors are defined 
by, 

C . — raia,2a s £0,3,0,20,! _ 

~ Jj-J J - -^-C( lj2 ,3)C( 3 ,2,l) , 

CyW = , r ^ c (l,2,3) c (3,4,5) c (5,6,7) c (7,8, l) c (9,2,10) c (10,4,ll) c (ll,6,12) c (12,8,9) , 

r 1 

^F( B ) — AT ^ C (1,4,3) C (3,2,5) C (5,6,7) C (7,8,1) C (9,2,10) C (10,4,11)C(H 1 6,12)C(12,8,9) > 
Cvi = CaC v (A) , 

Cy 2 = CaC V (b) , 
Cv 8 = CaC V (b) , 

d abcd d abcd = 1 NaCa (C vW + 2C V(B) ) , (Bl) 

where Na = Tr^l) is the number of gluon states, or the dimension of the adjoint represen- 
tation. The structure constants are written as 

c {i ,j,k) = f aiajak = Tr([T a \ T aj ] T ak ) , (B2) 

and i, j, k label internal or external lines. 

For G = SU(N C ), the quantities defined in eq. fIBip evaluate to, 



C A = 


2N C , 


Cy(A) = 


N C (N 2 + 12) , 


Cy(B) = 


12iV c , 


CVi = 


2N 2 (N 2 + 12) , 


CV 2 = 


2AN 2 C , 


Cv s = 


2AN 2 C , 


^abcd^abcd 


l(N 2 -l)N 2 c (N 2 + 3Q) 



(B3) 



Note that with our normalization (I5.3P of the generators T a for SU(N C ), namely Tr (T a T b ) 
5 ab , the value of Ca is twice the more conventional value. 
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For four-point amplitudes in any purely adjoint gauge theory, with gauge group G, we 
take the basis of color tensors through four loops (L < 4) to be {bf^}, where % is a label 
distinguishing color structures of same loop order. An explicit choice of basis is, 



bf 


— c (l,2,5) c (5,3,4) , 




bf 


= C( 2 ,3,5)C(5,4,1) , 




bf 


= C(i )5| 8)C(2,6,5)C(3,7,6)C(4,8,7) , 




bf 


= C (1,7,5) C (2,6,7) C (3,11,9) C (4,8,H) C (6,9,10) C (5,10,8) j 




6 (2) 


= C(2,7,5)C(3,6,7)C(4,11,9)C(1,8,11)C(6,9,10)C(5,10,8) , 




bf 


= c (l, 5,6) c (2,9,5) c (3,13,14) c (4,ll,13) c (6,7,8) c (8, 10,11) c (7,9,12) c (10,12,14) > 




bf 


— C (1,5,6) C (2,7,5) C (3,17,14)C(4,16,17) C (6,9,10) C (7,8,9) C (8,11,12) C (10,12,13) C (11,14,15) C (13,15,16) j 




bf 


= C(2,5,6)C(3,7,5) C (4,17,14)C(1 I 16,17) C (6,9,10)C(7,8,9) C (8,11,12)C(10,12,13) C (11,14,15) C (13,15,16) ■ 


(B4) 



b\ and b 2 are the s- and t-channel four-point trees, respectively; b\ is the one-loop box; 
bf are the s-channel ladder graphs; bf are the t-channel ladder graphs. 

We have demonstrated that this basis is complete for four-point amplitudes with only 
adjoint states and interactions proportional to f abc . We have done so by explicitly solving 
the Jacobi identities for all color tensors appearing at < L < 4 for a generic amplitude, and 
converting the reducible tensor structures to lower-loop color factors multiplied by group 
invariants. (Typical reducible tensor structures are color graphs with two- and three-point 
subgraphs.) Thus, any L-loop color tensor can be expressed in terms of a linear combination 
of the basis tensors bf , . . . , bf . Through four loops the pattern is that the number of new 
independent (irreducible) L-loop tensor structures are: two for even L, and one for odd L. 
The full basis then consists of the independent L-loop basis tensors, plus all the independent 
basis tensors for lower loop orders, which gives a total of (|L + 2) basis tensors for even loop 
order, and |(L + 1) for odd L. We believe that this pattern continues for L > 4 for the case 
of four-point tensors, but we have not proven it. 

In the standard case of G = SU (N c ) the basis color tensors are, 




N c { Tri234 + Tr 1432 ) + 2( Tr i2 Tr 34 + Tr 14 Tr 23 + Tr i3 Tr 24 ) , 



Tr 1234 + Tr 



Tr 1234 + Tr 



1432 



1432 ~ 



- Tr 



Tr 



1324 



1243 ~ 



- Tr 



Tr 



1342 , 



1423 > 
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bf ] = (N 2 + 2)( Tr 1234 + Tr 1432 ) + 2 Tr 1243 + 2 Tr 1342 - 4 Tr 1423 - 4 Tr 1324 + 6N C Tr 12 Tr 34 , 
b? = (N 2 C + 2) ( Tr 1234 + Tr 1432 ) + 2 Tr 1423 + 2 Tr 1324 - 4 Tr 1342 - 4 Tr 1243 + QN C Tr 14 Tr 23 , 
&f > = (AT3 + 2iV c )( Tr 1234 + Tr 1432 ) + 2iV c ( Tr 1243 + Tr 1342 ) + (14iV c 2 + 8) Tr 12 Tr 34 

+ 8 Tr M Tr 23 + 8 Tr i3 Tr 24 , 
bf ] = (iV c 4 + 2N 2 C + 8)( Tr 1234 + Tr 1432 ) + (2N 2 C + 8)( Tr 1243 + Tr 1342 ) 

- 16( Tr 1324 + Tr 1423 ) + (30iV c 3 + 24iV c ) Tr 12 Tr 34 , 

4 4) = (^c + 2iV c 2 + 8)( Tr 1234 + Tr 1432 ) + (2N 2 + 8)( Tr 1423 + Tr 1324 ) 

- 16( Tr 1243 + Tr 1342 ) + (3(W 6 3 + 24iV c ) Tr 14 Tr 23 . (B5) 



Using the Jacobi identity, the four-loop four-point color tensors Cj appearing in eq. (14. ip 
can be reduced to the eight-dimensional basis ( 1B4I) . For completeness, we first give the 
corresponding reductions of the one-loop four-point tensor, 



C 



'box 



(B6) 



the two-loop ones, 




(B7) 



and the three-loop ones, 



= &i 3) , 
cv> = b? - \C A bf\ 
c< c > = bf - \C A bf\ 
C® = b? - C A b? + \c%t}\ 
C w = -lb^ + p A b? + iC vm bf\ 

= -±bf + \c A bf - \C 2 A b? + \C v ^bf\ 
C« = -\b® + \C A b? + \C v ^bf\ 
C*( h ) = ±bf — j^Cj^b^ + ^C A b^ — \C V (B)bf\ 

= \C A b? - \C A bf. 



(B8) 
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The four-loop reduction gives, 



C 2 

c 3 
c, 
c 5 
c 6 
c 7 

c 9 

Cio 
di 

Cl2 
Cl3 
C*14 

C*15 
Cl6 
Cyj 

Cis 
Cig 

C*20 
C*21 
C22 
C23 
C24 
C*25 



('1 



6 



(4 



(i 



(4 



_Ia( 4 
2 U 1 

2°1 

l/,( 4 
2 



2°1 



2°1 
2°1 
2°1 



V 1,(3) 

ir» 1.(3) 

, 3 r l(3) , 1^ 7(0) 

+ \C A bf-\C\b^ + \C V2 bf\ 

, 1.(3) , 1^ l(0) 

+ Ca&i 3 — \c 2 A bf\ 

T 2°A^i — 8 A 1 ' 
(3) . 1^ 7(0) 



.I?/ 4 
2°1 



2 

3^ l(3) 



1l(4) l r 7(3) 1^ ,(0) 



bf - \C A bf + |Ci6^ 



l(3) 



y( B ) ( 



°1 + To ^^"l 7 



.(0) 



-&S 4) 


+ 


§crf> 


- 3 r 2 /) (2) 

8 U A°1 


r h {1) J- 5 r ^(°) 


-&! 4) 


+ 




1^2 7(2) 

2W1 




-&i 4) 


+ 




8°A°1 





2°1 



^i 4) - \c A bf 



\C A b^ + gC^fo^ + \Cy 2 bf\ 



(0) 



1^2,(2) 
4°A°1 



TnOy^ 



(0) 



-K» + 


\C A bf - 


1 r 2 b m - 


i r 1,(0) 
3^ 2 "l , 




-w + 


\C A bf - 


\c\bf - 


C V (B)b^ 


, 1 r 1.(0) 


+ 


\C A bf + 


\c\bf + 


I 8°A 




-W + 


±C A b? - 


\c\b? - 


Cy(B)b^ 


, I/O 7(0) 













TOy 2 1 



(0) 



1 r -L- 1/^2^(2) 1 ,(' 



(0) 
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^26 — I^AOi — ^A°l ~ 12°^2°1 > 

cw = |&i 4 — Icu&i 3 + §ci&i 2 — |cv 2 fri°\ 

^ - l^ 4 ) 1|,( 4 ) x V ^ 11^2,(2) , 5^2,(2) \ r ,(1) 7^ l(0) ,5^ iSP) 

n - U( 4 ) 11.(4) , 11^ l(3) 17^2^(2) 5^2^(2) i r ,(1) 13^ i(0) , 5^1 1,(0) 

^ — 1 5p 1,(3) 3ai2l(2) 1^ i(l) , 5 ^ iM \ r> 1,(0) 

^30 — — 2°1 s^^ ! _ 16 U A°1 ~ 2 L, l ,(B)fi i + 24°^2°1 — (j^VjA > 

^ _ U( 4 ) ^>r _l_ 13/^2 l(2) 5/t2i(2) , l-o l(1) ,19/^ 7,(0) 5^ l(0) 

^31 — 2°2 ~~ 8° A °1 + 16°A°1 ~~ 8°A°2 + 2°V( B >°1 + 24°^2°1 — Q L/ V 2 2 > 

^ - A( 4 ) 1 9^ 7,(3) , 35-^27,(2) 17^27,(2) ^ ,(1) 7 r 7,(0) 5>o l(0) 

^32 — 0\ + 2°2 ~~ I^AOi + 2i<^A l ~ 24 U A°2 + U V( B )°1 + 12°V 2 1 — 6 O ^2 2 5 

r> — U( 4 ) _L 5 r< A® 13^2^(2) 5^2^(2) 1^ t(1) z,(0) , 5-0 l(0) 

r _ 17,(4) , 7,(4) 9 r 7,(3) , 9r 2 h (2) 5 r 2 h (2), r ,(1) , 17 ^ l(0) 5^ 1,(0) 

^ — 1?,( 4 ) _u 1,(3) 35^27,(2) 13^27,(2) l^y ,(1) 5 r 7,(0) , 5 r 7,(0) 

^ _ 17,(4) 1 r 7,(3) , 15^2 1(2) 5 r 2 h {2) , 1 ^ ,(1) 7,(0) l(0) 

^ — l r _L 7^2i(2) l r 2i(2) l^y 7(1) 1 i r 7(0) 

C38 = + 4^ — 2Ca&i 3 ^ + H^a^i 2 ^ — I^i4 2 ^ + Cy^b^ + |Cy 2 fci°^ — |Cv 2 6 2 °\ 

r - J- ^( 4 ) /, (3 ^ J_ 17^7,(2) 5 r 2r(2) l r ,(1) l 7r y 7,(0) 5^ 1,(0) 

_ 17,(4) 7,( 3) -i- 1 r 2 h^ _u ^ ?,(°) 
C41 = — |C J 46p' ) + \C\b^ + Cy(B)6i 1 ' ) — ICva^i ^, 

r - 17,(4) 1 3 r 7,(3) 1^27,(2) , 1 7,(0) 
O42 — — 2°1 + I^AOi — 4^A°1 + j2 Uv 2°l ' 

^ -17,(4) 3 ^ 7,(3) , l r 2 h (2) , r 7(1) l r ,(0) 

— ir /,(°) 

O44 — -gO^Oj , 

r< -17,(4) 3 r 7,(3) , 1^27,(2) , ^ 7,(1) i r 7,(0) 
C45 — 2°1 — 4°^°1 + 4W1 + °y( B )°l — 4 L 'V 2 °1 ) 

C 46 = — I&t 4 ^ — 6 2 4 ^ + |Ca^i 3 ^ — 2C A bf^ + |C^6 2 2 ' ) — C V (n)b^ — j|Cy 2 6^ + |Cy 2 6 2 °\ 

C47 = |&^ 4 + 6 2 4 ^ — fCA^f' 1 + 2C A 6[ 2 ' ) — \C 2 A b^ + CytB)^! 1 ^ + C'v 2 &i ' > — §CV 2 &2°\ 
r - i/,(4) , 3 ^ 7,(3) \ r 2 h {2),\ r ,(0) 

O48 — — 2°1 + 4°^ 1 — i^A^l + l2 O V 2 l 7 

C49 = l/j^ + i/jg 4 '' — |C J 4&i 3 ' ) + fCl^l 2 '' — |C A /J 2 2 ' ) + Cy(B)b^ + |Cy 2 6f' ) — fCva^, 

C 50 = 0. (B9) 

In these equations, Cy 2 is the Casimir for the four-loop vacuum graph V 2 and C V (b) is the 
one for three-loop vacuum graph V^ B \ as defined in eq. (IBlj) . 

Finally, because the color tensors may appear with arbitrary permutations of the 
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external legs, we also need to give the reductions of the basis elements for arbitrary permu- 
tations. In the following, we quote particular permutations, as well as symmetries of the 
tensors that allow us to reduce all other permutations: 



S 0) (l,2,3,4) 


_ 6 (o) 


S 0) (2,3,4,l) 


_ 6 (o) 

— U 2 5 


S 0) (l,3,4,2) 


- b i0) - 


if' (a, 6, c, d) 


= -&f» 



do) 



(BIO) 



6^(1,2,3,4) = 6^, 
6^(1,3,4,2) = b^-\C A bf\ 

b^\a,b,c,d) =b?\b,c,d,a) = b?(a,d,c,b), (Bll) 



6^(1,2,3,4) 










6^(2,3,4,1) 


— U 2 5 








6S 2) (1,2,4,3) 


= &i 2) - 




• 




6^(2,3,1,4) 


= 4 2) - 


l/o2i(0) 
4 U A 2 


• 




6! 2) (1,3,4,2) 


= -b? 




^C A b^ — 


l r 2A0) 


6! 2) (1,3,2,4) 


= -b? 




fC^ 1 — 


1/021(0) 
4°A°2 


b^ (a, 6, c, d) 


= b?(b, 


a, d, c) = 


= b?\c,d, 


a, 6), 



(B12) 



6i 3) (l,2,3,4) 


= 6f 


5 




6^(1,2,4,3) 


= ff> 


1^3 1(0) 




^ 3) (2,3,4,1) 


= 6f 


-lC A {b?-b?)- 


W( B H°1 — °2 J) 


6^(2,3,1,4) 


= 6? 


-icM ] -bf)- 


<->y( B H°l - 2 J - gt>A°2 ) 


6^(1,3,2,4) 


= if 


-lC A (2b? + b?)- 




6S 3) (1,3,4,2) 


= lf> 


- lC A {2bf +6?))- 


VlClb^+Cy^bf -\C\bf 


b^\a, 6, c, d) 


= if » 


(6, a, d, c) = 6| (c, d, 


a, 6), 



(B13) 
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6S 4) (1,2,3,4) = 6 (4) 



1 > 



6S 4) (2,3,4,1) = 4 4 \ 
&< 4) (1,2,4,3) = -±C\b?\ 
6} 4) (2,3,l,4) = 6? ) -iCl6 3 W 

6S 4) (1,3,4,2) = -^ 4) - 4 4) + f«S 3) - ^C\bf + |(15CI - 16C v( b))6« 

— ^(SCl + 56(7^)6^ + yCV 2 & 2 ° , 
6S 4) (1,3,2,4) = - bf + \C A b® - 2 -{C\bf + |(15CI - 16^)^ 

-i(3Cl-160Cy 2 )fef -fcrf, 
6i 4) (a,^c,d) = &S 4) (6,a,cZ,c) = bf{c, d, a, b). (B14) 
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Appendix C: Numerator and Color Factors of the parent integrals 

In this appendix we collect all 50 numerator factors of the parent integrals, displayed in 
figs. I2D1I251 These polynomials are, 

7V 6 , N 7 , iVg, iV 9 , N 10 , N u = s? 2 s 45 , 
N12 = s1 2 ssq , 
N13 = - s\ 2 s m , 

N 15 , N 16 , N 17 = - s 12 s 2 56 , 

s 12 s 35 s 46 , 
-^19 = — S12S56S47 , 
N 2 Q, N 2 1, ^22, N 23 = Si2S 56 S 4 7 , 
N24 = 512(335338 — S37S36) , 

N 2 5 = S 12 (s 45 S 3S - S12S37) - ^Sl2S45 ~~ ^8 S 12 , 

N 2 6 = sl 2 S 8 9 ~ S 12 S 35 S 67 + /9S12T35 + ^S1 2 (S4,10 ~ ^ll) ) 

N 27 = Si2(Si2S78 + £4559,10) — S12S45 — 512(512/11 + S45/12) ) 

N 2 8 = S12S17S36 + S23S15S38 — -/9S12S23 — ^8 S 12 ~ ^7 S 23 , 
-^29 — S12S17S68 + S23 s 38 s 15 — ZgSi2S23 + ^8 S 12 ~ ^7/9^23 , 
^30 = S23S15T29 — s 13 s 26 r 19 — s 12 s 15 s 26 ~ ^)(^5 S 13 + ^l s 23) 5 

N31 = - S12S36S56 - S12S26S68 - S12S36S25 - S23S38S25 + /e s i2( s 57 + s 5 8 + s 8 g) - lts X2 {ll + If) , 

^32 = S 56 (si2S5 6 - S23S35 - llsi 2 ) , 

^33 = S 12 (si 5 S 7S - S29S56) + S 2 3S56(-S9,10 ~ ^lo) + /s( S 12 s 29 ~ /?2 S 23) , 

^34 = Si 5 (si3S 56 - S12S78 - ljsi 3 ) - Si 6 (s23-S 5 6 ~~ S 12-S79 - + ^Sl2Tl,lo , 



96 



^35 = S 2 3(si0,ll( s 17 - ll) ~ + Si 2 Si 5 (s 2 3 ~ 5 7 8 + Zg) 

+ Sl 6 (si2S 7 9 + Z10S23) - ^(Sl 2 Sl,10 - ^12 S 12 + ^10 S 23) , 

^36 = Si2S 4 7(s 2 3 - -§67 + ^0) + S^^^IO - S23S47S89 + lnS 12 (s e7 - S23) 
+ s 12(^7 _ ^6,10 — ^7 S 15) + ^11^12 s 23 + 2^5 S 12 S 23 > 

^37 = S 4 6(si 2 S28 - S2 3 S 67 ) + S 45 (si 3 S 37 - S 2 3S 58 ) + S23S45S46 + igSi 2 Sl3 ~ «23^6 > 

^38 = S56(si2S78 + S23S56 — S12S23 — Sj^Slo.ll ~ s 23^g) + S 12 S 10,11 + S 13(^12^14 + ^13^15) 1 

-^39 — s 12('Sl2S6 ) T4 — s 29 s 38 ~~ s 35 s 16,17 + S 1,15 S 4,17 — s 2,12 s 15,16) + s 13 s 15,16 s 16,17 

+ Sl2(^i4 s 5,T0 + ^13 S 7,T2 _ ^9 S 8,14 _ ^8 S 9,T3~) ~ S 13^5 S 12,14 + ^12 S 5,13) 

_L e r/ 2 / 2 — l 2 ! 2 \ -L q |'/ 2 7 2 _l/2;2 _ /2/2 , ;2/2 1/2/2 _ /2 /2 _l_ ] 2 ] 2 \ 
~r ^23^l5 l 17 5 12/ ' 12\ 8 9 7 10 l 7 l 13 ~r ^'lS "r l 9 l 14 l 10 l 14 1" L l2, l lA) 1 

N m = - s 2 12 s 5G , 
N41 = S12S45 , 

^43 = - S12S35S46 , 

N44 = 512(545557 — 53555s) , 

^45 = ^12(545578 + 5355g 5 io — 545S35 — l\iS^ — ^545) j 

^46 = 5 68 (5 2 3545 - 5i 2 5 46 ) + S7 8 (5i25 4 7 ~ 5i 3 5 4 5) + ^6 S 13 ~ I7S23) + ^7 S 13 ~ ^23) , 
^47 = Si3S6,lo(568 — s 12 ~ ^11) ~ 52358,10(568 ~~ s 12 ~ ^11) 
+ 5i2(s 7 95 6 ,10 - S595 8 ,lo) + 2s 12^6 _ l l)( l 12 ~ l lo) > 

N48 = 5i 2 (s 2 ,10539 ~ S47S18 + S 2 ,1 S 5 9 + S 3 9S6,10 + S23S6.ll) ~ 523^57568 ~ 513^5956,10 
+ ^6( S 12 S 35 + S1254J2 ~~ s 23 s 59) + ^s(Sl2526 + S 12 S 1,TT — s 23 s 6,lo) 
+ ig(si2S 12 j3 - Si 3 Sio,ll) + ^ ( s 12Sll,n ~ S 13-S9,12) 
- /i35i25 n) i4 - /i4 s 12Si2,T3 + ( s 13 ~ 2sl 2 )ZgZ? 
+ ^23(^5^6 ~~ ^8 + ^7 + W'l) + s 12^14 - 512^5^6 

+ Sl2(-^5^8 + ^9 ~ Mil ~~ ^15 ~ ^9^15) 

_L_ o f_/ 2 /2 , /2/2 _ /2/2 _ /2/2 /2 /2 \ 

T 5i2l, '6 l 7 "r ^'lO (, 6 t 12 l 6 l 16 l 10 l 16^ 

+ s 23(^9^12 + ^10^11 ~~ ^9 ~~ ^lo) + s 13(^ll + ^12) ; 
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^49 = Si 2 (s47S 5il 2 - Si 9 S 36 - S48S36) + 523(548^6,11 ~ Si 5 S 3 ,10 ~ S15S47) 
+ ^( s 23 s 7,12 — S23 s 4,15 — s 13 s 10,ll) + ^6 ( S 12 s 8,ll ~~ s 12 s 4,l5 ~~ s 13 s 9,12) 
+ ^(523^3,15 - S12S38 + S 2 3S 6 ,lo) + ^lo( S 12Sl,15 - S 23-Si7 + S12S59) 
+ ^13( S 12S23 + S12S38 — S23 s 6,ll) + ^i 4 (s23 s 12 + S23 s 17 — 512^5^2) 
+ ^ll S 23( s 4,12 — s 6,lo) + ^12 S 12(S4,11 — S59) 



Sl2S23Sll,12 



_L /;2;2 , ,2,2, ,2,2 , ,2 ,2 , ,2 ,2 , t2/2 _/2,2 /2 ,2 \ 

T &13^7<'8 t 5 ,, 8 L 6 L 7 ' l ll l 12 ~r none L 9 L 17 t 9''12 nCTllJ 

+ ^H^O + mi + lf 3 ~ &) + M + i§ - i? + + mi + 4)) 

+ S23(-^9 + ^(^13 + ^14 ~~ ^9) + ^ll(^10 + l\ ~ l\ + ^13) + ^(^10 + ^17)) 



_L_ e (I 2 /2 _ /2/2 _ ,2 ,2 _ /2 /2 ,2 ,2 \ . //2 ,2 _ ,2/2 _ ,2/2 _ ,2/2 _ ,2 ,2 \ 

T ^12^i2 t 13 l 8n3 l 10 l 13 l 10 l 14 t 13 t 17/' "r i, 23^ll t 14 l 7 l 14 f 9 l 14 f 9 l 13 'uneJ ' 

— s 12'528 s 4,12 — Sl2S37 s l,ll — S23 s 16 s 3,10 + S23 s 25 s 49 + ~ S12S23 ($13,15 — s 13,14) 



We note that graph 50 has a vanishing color factor and therefore does not contribute; the 
integral also vanishes under integration. 

We now give the tensor color factors for the 50 four-loop integrals describing the ampli- 
tude. Although it is simple to read off these color factors from the graphs, we include them 
here because they contain signs. The color-factor sign (and hence the sign of the numerator 
factor multiplying it) depends on how the graph is drawn. The signs here are properly 
correlated with the signs of the numerator factors given above, as well as with the graphs 
in figs. EUH251 We have, 

— C (1,5,6) C (2,7,5) C (3,17,14)C(4,16,17) C (6,9,10) C (7,8,9)C(8,11,12)C(10,12,13) C (11,14,15) C (13,15,16) j 

— C (1,5,6) C (2,7,5) C (3,15,16) C (4,17,14)C(6,9,10) C (7,8,9) C (8,11,12)C(10,12,13) C (11,14,15) C (13,16,17) > 

— C (1,5,6) C (2,9,8) C (3,14,15) C (4,17,16) C (5,7,8) C (6,9,10) C (7,11,12)C(10,12,13) C (11,16,14)C(13,15 ) 17) , 

— C (1,5,6) C (2,7,5) C (3,13,16) C (4,16,17) C (6,9,10) C (7,8,9) C (8,11,12)C(10,14,15) C (11,15,17) C (12,13,14) 5 
C5 = C(i i 5 j 6)C(2,9,8)C(3,13,16)C(4,16,17)C(5,7,8)C(6,9,10)C(7,ll,12)C(l0,14,15)C(ll,15,17)C(12,13,14) , 

— c (l,6,7) c (2,8,6) c (3,16,ll) c (4,15,17) c (5, 14,15) C (7,10,5) C (8,9,10) C (9,11,12)C(12,13,14)C(13 1 16,17) 1 

— C (1,6,7) C (2,8,6) C (3,15,13) C (4,12,17) C (5,11,12)C(7,10,5) C (8,9,10) C (9,16,15) C (11,13,14)C(14,16,17) 1 

— C(l,6,7)C(2,8,6)C(3,13,16)C(4,15,17)C(5, 14,15) C(7,10,5)C(8,9,10)C(9,11,12)C(H,17 I 16)C(12,13,14) , 

— C (1,6,7) C (2,10,9) C (3,17,11) C (4,15,16) C (5,14,15) C (6,8,9) C (7,10,5) C (8,11,12)C(12,13,14)C(13,17 I 16) , 



+ 312(^10 




(CI) 
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Cio 




c 


1,6,7) C (2,8,10) C (3,12,13) C (4,16,17) C (5,15,16) C (6,9,10) C (7,8,5) C (9,11,12) C (11,17,14)C(13 I 14,15) 5 


C u 


= 


c 


1,6,7) C (2,10,9) C (3,13,16) C (4,15,17) C (5,14,15) C (6,8,9) C (7,10,5) C (8,11,12)C(H I 17 I 16)C(12,13,14) 9 


Cl2 


= 


c 


1,7,5) C (2,9,7)C(3 I 15 I 16)C(4 I 6,15) C (5,8,13) C (9,11,10) C (10,12,8) C (11,16,17) C (12,17,14) C (13,14,6) 9 


C*13 


= 


c 


1,7,5)C(2,8,7)C(3,17,16)C(4,6,17)C(5,13,H) C (8,9,10)C(9,11,12)C(10,15,14)C(12,16,15)C(13,14,6) 9 


C14 


= 


c 


1,8,9) C (2,5,8) C (3,14,16) C (4,10,13) C (7,17,15) C (6,15,14) C (5,16,17) C (9,11,10) C (11,7,12)C(12,6,13) > 


Cl5 


= 


c 


1,7,8) c (2,9,7) c (3,15,17) c (4,14,16) c (8,12, 13) c (9,10,ll) c (10,6,15) c (ll, 17,5) c (12,5, 16) c (13,14,6) 9 


Cl6 


= 


c 


l,7,8)C(2,9,7)C(3,17,16)C(4,12,13)C(8,13,14)C(9,lo,ll)C(io,6,17)C(ll,16,5)C(12,5,15)C(14,15,6) 9 


C17 


= 


c 


l,7,8) c (2,8,9) c (3,14,16) c (4,ll,12) c (5,15,17) c (6,14,15) c (7,17,16) c (9,10,ll) c (10,6,13) c (12,13,5) 9 


Cl8 


= 


c 


1,8,9) C (2,5,8) C (3,6,12)C(4,16,14) C (7,17,16) C (13,15,17) C (5,12,11) C (9,10,7) C (10,11,13) C (6,14,15) 9 


Cl9 


= 


c 


1,8, 7) c (2,9,8) c (3,15, 11) c (4,16, 14) ^5,11,12) C(7 > i6,17) c (9, 10,5) c (10, 17,6) c (12, 13,14) C(13,15,6) 9 


C20 


= 


c 


1,8,9) c (2, 7,8) c (3,15, 13) c (4,ll, 12) ^7,10,11)^9,6, 17) c (10, 17,5) c (12, 13, 14) C(14,16,6) C (15, 5, 16) > 


C21 


= 


c 


1,8, 7) C(2,9,8)C(3,13,14)C(4,17,11)C(5, 14,15) C(6,12, 13) C(7,16,17)C(9, 10,5) C(10,11,6)C(12, 16,15) 9 


C22 


= 


c 


1,8,7) C (2,9,8) C (3,16,14)C(4,12,11) C (5,15,17) C (6,17,16) C (7,13,12) C (9,10,6) C (10,11,5) C (13,14,15) 9 


C23 


= 


c 


1,8,7) C (2,9,8) C (3,15,17) C (4,12,11) C (5,17,16) C (6,16,14)C(7,13 1 12)C(9 1 10,6) C (10,11,5) C (13,14,15) 9 


C*24 


= 


c 


1,10,7)C(2,5,10) C (3,16,14)C(4,13,9)C(5,11,6)C(6,17,16)C(7,8,12)C(8,14,15)C(11,12,13)C(15,17,9) 9 


C25 


= 


c 


l,10,6) c (2,5,10) c (3,14,16) c (4,17,8) c (5,ll,9) c (9,12,13) c (6,7,15) c (ll,15,17) c (12,8,16) c (13,14,7) 9 


C*26 


= 


c 


1,12,13) C (2,5,12)C(3,16,11) C (4,15,9) C (5,14,6) C (6,17,16) C (7,9,17) C (10,11,15) C (13,10,8) C (14,8,7) 9 


C27 


= 


c 


1,13,5) C (2,6,13) C (3,16,15) C (4,17,7) C (5,9,8) C (6,14,10) C (10,15,11) C (11,17,9) C (14,8,12)C(16,12,7) 9 


C*28 


= 


c 


l,10,ll) c (2,6,5) c (3,15,13) c (4,7,8) c (6,13,14) c (10,5,12) c (ll,17,7) c (12,9,17) c (14,16,9) c (15,8,16) 9 


C*29 


= 


c 


1,10,11)C(2,9,5) C (3,17,15)C(4,7,8) C (9,13,14)C(10,5,12)C(H,16,7) C (12,6,16)C(13,8,15)C(14,17,6) 9 


C30 


= 


c 


1, 10,11) c (2, 15, 13) C (3,6,5) C (4,8,7) C (6,13,14)C(10,5, 12) C(ll,17,9) c (12,7,17) c (14,16,8) c (15,9,16) 9 


Ga 


= 


c 


1,9,5) c (2, 16,10) C (3,17,14)C(4,7,8) C (5, 10,11) c (9,12, 13) c (ll, 15, 7) c (12,8, 14) ^13,17,6) C(i5,i6,6) 9 


C32 


= 


c 


1, 7,6) C(2,9,15)C(3, 17,16) C(4,5,8)C(5, 10,11) C(6, 12, 10) C(7,8,9)C(11, 14,17) C(12, 13, 14)C(13,15,16) , 


C33 


= 


c 


1 1 5 1 71 C(o 1 9 1 1 )C(3 7 q\C(A 1 n 1 3) Cf6 1 fi 1 41 Cfa q 1 n\ C(-\ n « filCn 3 <^ 71 Cfi 4 1 ^Cfi f! 1 1 1 71 , 


C34 




c 


1,16,14)C(2,13,7) C (3,11,9) C (4,8,12)C(5,14,15) C (6,17,16) C (10,5,8) C (10,9,6) C (11,12,13) C (15,17,7) 9 


C35 




c 


1,16,14)C(2,13,7) C (3,8,11) C (4,10,12)C(5,14,15) C (6,17,16) C (9,11,13) C (10,9,6) C (12,5,8) C (15,17,7) 9 


C36 




c 


l,13,9) c (2,8,6) c (3, 12, 7) c (4,15, 17) c (5, 14,15) C (8,10,12)C(9 ) 5,10)C(H ) 16, 14) ^13,6,11) C(i6,7,17) 9 
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C*37 = C(i i io,ll)C(2j l l2)C(3 i i3 ) 8)C(4 i i6,14)C(6,17,16) c (7,8,9) c (9,15,17) c (10,12,6) c (ll,5,13) c (14,15,5) > 
C38 — C (1,16,12) C (2,13,8) C (3,14,7) C (4,17,15) C (5,15,10) C (6,10,14) C (9,7,17) C (12,5,H) C (13,11,6) C (16,8,9) , 
C39 = C (1,5,9) C (2,13,10) C (3,14,7) C (4,12,8) C (5,15,6) C (9,10,17) C (13,6,16) C (15,8,7) C (16,14,11) C (17,11,12) ) 
C40 = C(i i 7 i 8)C(2,9,7)C(3,14,15) c (4,16,14) c (5, 17,16) C(8,12,13) c (9,6, 17) C (10,5,12) C (10,6,11) C (H, 13,15) ; 
C41 = C(i i 6 j 7)C(2,5 l 6)C(3 i i6,15) c (4,10,8) c (5,8,9) c (7,15,13) c (9,14,17) c (10,ll,12)C(ll,l3,l4)C(i2,16,17) 1 
C42 = 0(1^1^)0(2,5,11)0(3^6,12)^4,9,7)0(5^3^4) 0(7,12,8)0(8, 15, 13) c (9,6,10) c (10,15,17) c (14,16,17) > 
C*43 = 0(1,7,8) 0(2,5,7)0(3,15,9)0(4,16,10) 0(5,6,9)0(6,10,11)0(8,12,17)0(11,12,13)0(13,14,15)0(14,16,17) 9 
C44 = C(i,6,l0)C(2,10,5)C(3,7,13)O(4,8,12)C(5,l6,14)C(6,ll,9)O(8,14,15)O(9,17,15)O(ll,12,13)O(16,7,17) , 
C45 = 0(1,13,6)0(2,5,13)0(3,9,14) C(4,i6,8)O(7,6,10)O(7,ll, 14)0(10,12,16)0(11, 8, 15) 0(12,9,17)0(15,5,17) , 
C*46 = C(i, n,6) 0(2,7,12) 0(3,5,9) 0(4,10,13) 0(5,13,8) C(6,16,14) 0(8,16,17) c (10,14,15) 0(11,9,12)0(15,7,17) , 
C47 = C(i,6,5) 0(2,8, 7)0(3,12,13)0(4,9, 10) 0(5,7,11)0(6, 14,15) 0(11, 13,9)0(14,10,16) 0(15, 12, 17) 0(16,8,17) > 
C48 = C(i,9,i5)C(2,5,13)C(3,14,6) 0(4,10,16)0(5,6,17) 0(7,12,10)0(8,14,16)0(9, 11,8)0(15, 7,13)0(17,12,11) , 
C49 — C(i,i6,ll)C(2,6,5)C(3,12,17)C(4,10,9)O(6,17,8)C(7,9,14)O(8,13,10)O(ll,13,15)O(14,12,15)O(16,5,7) , 
C50 = C(i,i2,5)C(2,7,6)C(3,9,8)C(4,10,ll)C(5,13,6)C(7,8,14)C(l2,15,ll)O(13,16,17)O(14,16,15)O(17,10,9) • (C2) 

Group indices ai, a-i, as, a\ correspond to external color, and are not summed over. Internal 
indices 05 through a i7 are to be summed over. Thus the suppressed tensor structure of each 
color factor is Cj = Cf 10,20,3 "" 1 . It turns out that C 50 vanishes identically for all gauge groups 
(because of constraints from the Jacobi identity). 
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